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ABSTRACT 

We present time-resolved broad-band observations of the quasar 3C 279 obtained from multi- 
wavelength campaigns conducted during the first two years of the Fermi Gamma-ray Space Telescope 
mission. While investigating the previously reported 7-ray/optical flare accompanied by a change in 
optical polarization, we found that the optical emission appears delayed with respect to the 7-ray 
emission by about 10 days. X-ray observations reveal a pair of 'isolated' flares separated by ~ 90 
days, with only weak 7-ray/optical counterparts. The spectral structure measured by Spitzer reveals 
a synchrotron component peaking in the mid-infrared band with a sharp break at the far-infrared band 
during the 7-ray flare, while the peak appears in the mm/sub-mm band in the low state. Selected 
spectral energy distributions are fitted with leptonic models including Comptonization of external 
radiation produced in a dusty torus or the broad-line region. Adopting the interpretation of the po- 
larization swing involving propagation of the emitting region along a curved trajectory, we can explain 
the evolution of the broad-band spectra during the 7-ray flaring event by a shift of its location from 
~ 1 pc to ~ 4 pc from the central black hole. On the other hand, if the 7-ray flare is generated instead 
at sub-pc distance from the central black hole, the far-infrared break can be explained by synchrotron 
self-absorption. We also model the low spectral state, dominated by the mm/sub-mm peaking syn- 
chrotron component, and suggest that the corresponding inverse-Compton component explains the 
steady X-ray emission. 

Keywords: radiation mechanisms: non-thermal — galaxies: active — galaxies: jets — quasars: indi- 
vidual (3C 279) — gamma rays: galaxies — X-rays: galaxies 
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1. INTRODUCTION 

Blazars are Active Galactic Nuclei (AGN) character- 
ized by highly luminous and rapidly variable continuum 
emission at all observed bands. The most commonly ac- 
cepted scenario has their broad-band emission Doppler- 
boosted by a relativistic jet p ointing close to our line of 
sight (e.g.. lUlrich et al.l 119971 ). While the jet emission 
usually dominates the observed broad-band spectrum, 
the optical/ultra- violet (UV) and infrared (IR) spectra 
often also reveal signatures of the central engine: broad 
emission lines, and in some cases, quasi-thermal opti- 
cal/UV emission and IR dust emission, indicating the 
presence of an accreting supermassive black hole. Most 
viable current models for the origin of such jets involve 
conversion of the gravitational energy of matter flowing 
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onto the black hole to the kinetic energy of the relativistic 
outflow or tapping the rotation energy of a spinning black 
hole. However, the conversion process itself is not well 
understood, and many additional questions regarding the 
dissipation region of the jet's energy into radiation and, 
in particular, its location remain unanswered. 

Major advances in understanding of blazars came as 
a result of the discovery by the EGRET instrument on 
board Compton Gamma-Ray Observatory (C'GRO) that 
they are strong 7-ra y emitters, with 7-r ays dominat- 
ing radiative output (jHartman et a.1.1 Tl999T ) . With this, 
multi-band observations including the 7-ray band hold 
the promise of answering many outstanding questions re- 
garding the structure of t he relativistic jets of blazars. 

3C 279 (z = 0.536; ILvnds et all [19651) is in fact 
one o f the first 7-ray blaza rs discovered by EGRET in 
1991 (Ha rtman et al.|[j_992T) . The 7-ray signal had been 
significantly detected in e ach observation by EG RET 
since its discovery (see, e.g.. lHartman et aT1l200 1a). with 
the flux having ranged over roughly 2 orders of magni- 
tude, from ~ 10~ 7 up to ~ 10~ 5 photons cm" 2 s -1 above 
100 MeV (jMaraschi et al.lH99l IWehrle et alj|l998l ). and 
a factor of 2 variation on timescales as short as 8 hrs. 
The photon in dex in the EGRET 7-ra y band ranged 
from 1.8 to 2.3 (jNandikotkur et all 12001 . On a few oc- 
casions, 3C 279 was also detecte d at lower energies by 
CGR O's OSSE (50 keV-1 MeV) (IMcNaron-Brown et al 
1995) and CQMPTEL (0 .75-30 MeV) (|Hermsen et al 



19931 ICollmar et al.ll200lT ) instruments, indicating that 



the 7-ray emission forms a broad peak in the vF v rep- 
resentation. In 2008 July, the AGILE satellite ob- 
served a 7-ray flar e associated with th e source with 
11.1 a significance (|Giuliani et al.l [20091. with an av- 
erage flux above 100 MeV of (21.0 ± 3.8) x 10" 7 
photons cm~ 2 s _1 and the photon index of 2.22±0.23 be- 
tween 100 and 1000 MeV. In the very-high-energy (VHE) 
7-ray regime above 100 GeV, the imaging atmospheric 
Cherenkov teles cope MAGIC dete cted flares twice in 
2006 February (lAlbert et alj I2009D and in 2007 Jan- 
uary (jAleksic et al.l l2011h . which made this source the 
most distant currently known VHE 7-ray emitter. 

Optical and UV observations of the source in relatively 
low states - when the jet emission was relatively faint - 
allowed a study of the accreting black hole and the associ- 
ated accretion disk. The luminosity of the ac cretion disk 
was e stimated to be Lu ~ 2 x 10 45 erg s _1 (jPian et al.1 
1999). The mass of central supermassive black hole was 
estimated to be in the range of ( 3-8) x lQ^M^ using the 
luminosity of optica l broad line (|Woo &: Urrvl [2002t ) or 
the Hp line width (|Gu et al.1 120011 ) . Those values are 
similar to th e estimates based o n the luminosity of the 
host galaxy (jNilsson et al.ll2009[ ). 

3C 279 contains a compact radio core, associated with 
time-variable jet-like structure. Radio observations at 
43 GHz by Very Long Baseline Array (VLBA) between 
1998 March and 2001 April revealed superluminal mo- 
tion of the jet with apparent spee ds that range from 5c 
to 17c porstad et all 12004 I2005D . Those observations 
also allowed an estimate of the Lorentz factor of the jet 
flow of Fj = 15.5 ± 2.5 and of the viewing angle of the jet 
Oo = 2.1 ± 1.1 deg, corresponding to a Doppler beaming 
factor of S — 24.1 ± 6.5. A change of the trajectory of 
a jet component has also been reported in radio obser- 
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vations with VLBA (|Homan et al .1120031) . Those authors 
estimated the jet component to be moving with a Lorentz 
factor Tj > 15 at an initial viewing angle of < 1°. 

The broad-band spectral energy distribution (SED) 
of the source is characterized by a two-bump struc- 
ture, similar to many other 7-ray blazars. In the 
context of widely accepted leptonic models, the lower- 
frequency bump, peaking at the far-IR and extend- 
ing to the extreme UV band, is commonly ascribed to 
synchrotron radiation from relativistic electrons in the 
jet. The second bump, spanning from the X-ray to 
the 7-ray band with a peak in the MeV-GeV range, 
is believed to be generated via inverse-Compton scat- 
tering, presumably by the same population of parti- 
cles that radiate at lower energies via the synchrotron 
process. The seed photons for the Compton scattering 
can b e synchrotron ph otons (synchrotron self-Compton : 
SSC, IMaraschi etafl Il992t iBloom fc Marscherl 119961) . 
accretion disk photons (external Compton sca t tering 
of direct disk radiation: ECD, iDermer et al.l 119921 : 
iDermer fc Schlickeiserlll993l ) and accretion disk photons 
re-scattered by the broad-line region clouds /intercloud 
medi um (ECC. ISikora et al1ll994HBlandford fc Levinsonl 
1995), or infrared radiation from a torus located be - 
yond the broad- line region (ERC-IR, [Sikora et al"1ll994l ). 
Specific to 3C 279, multiwavelength snapshot observa- 
ti ons for several epochs in cluding 7-rays were presented 
in IHartman et al.l (|2001af ) . Those authors explained the 
overall spectra using the leptonic model, where the X-ray 
photons are mainly produced by SSC, and both ECD and 
ECC contribute to the 7-ray emission. Spectral variabil- 
ity was explained by variations of the bulk Lorentz factor 
of the jet, accompanied by changes in the spectral shape 
of the electron distribution. 

The optical variability of 3C 279 is ext reme: in 1937, it 
show ed optical (B) magnitude of 11.27 IjEachus fc Lillerl 
I1975T ). making it one of the most luminous active galax- 
ies ever recorded. The strong variability recorded in all 
bands provides an opportunity to establish the relation- 
ship between emission in those bands, and thus can be 
used to constrain theoretical models of physical regions of 
the jet responsible for such emission. Many such multi- 
wavelength campaig ns have been condu cted (see, e.g., 
Maraschi et alJll9M iWehrle et aflll998t iLarionov et al.l 
2008t iCollmar et alj |2010| ) but thev have not revealed 
a simple relationship between the variability in various 
bands: radiation in different spectral regimes does not 
always rise and fall simultaneously, although the periods 
of increased rapid activity in all bands seem to last for 
several months, and take place when the source i s rela - 
tively bright. A recent paper bv lChatteriee et al.l ([20081 ) 
presents the results of the monitoring of 3C 279 for 11 
years in radio, optical and X-rays, and discusses the de- 
tails of the jet structure based on multi-band correlation 
studies. However, due to the lack of deployed instru- 
ments, long term monitoring observations could not in- 
clude the 7-ray regime, where the source often shows 
stronger variability than in other bands. 

The launch of the Fermi Gamma-ray Space Telescope 
on 2008 June 11 has rejuvenated multi-b and studies of 
blaza rs. The Large Area Telescope (LAT; lAtwood et al.l 
2009) instrument on Fermi can monitor all 7-ray sources 
on the sky with its a wide filed of view and a much larger 
effective area compared to earlier 7-ray missions. Taking 



advantage of this new instrument for 7-ray observations, 
we have organized intensive multiwavelength campaigns 
for 3C 279 from radio to the high-energy 7-ray energy 
ranges. Many ground-based telescopes (cm, mm, near- 
IR and optical) and various satellites (IR, UV, X-ray, 
hard X-ray and high-energy 7-ray) participated in this 
ca mpaign. We reporte d the first results of the campaign 
in lAbdo et alj (|2010al hereafter Paper I), where we dis- 
covered the dramatic change of the optical polarization 
coincident with the 7-ray flare. Here, we provide details 
of the multi-band observations and the interpretation of 
those data for the 2-year interval between 2008 August 
and 2010 August. In Section [21 we present and briefly 
discuss the features of the LAT 7-ray data; in Section 
02 we present the data in lower energy bands. Section [4] 
highlights the features of time series measured in various 
bands including their cross correlations, and the general 
properties of the broad-band spectral energy distribu- 
tion. In Section we provide viable emission models for 
the source in the context of leptonic scenarios. 

2. Fermi-L AT DATA AND RESULTS 

Ferrm-LAT is a pair-production telescope with large 
effective area (8,000 cm 2 on axis at 1 GeV for the event 
class considered here), and large field of view (2.4 sr at 
1 GeV) , sensitive to 7 rays in the energy range from 
20MeV to > 300 GeV. Informatio n regarding on-orbi t 
calibration procedures is given in lAbdo et alj ([2009a). 
Fermi-LAT normally operates in a scanning 'sky-survey' 
mode, which provides a full-sky coverage every two orbits 
(3 hours) . For operational reasons, the standard rocking 
angle (defined as the angle between the zenith and the 
center of the LAT field of view) for survey mode was 
increased from 35° to 50° on 2009 September 3. 

2.1. Observation and data reductions 

The data used here comprise 2-year observations ob- 
tained between 2008 August 4 and 2010 August 6 (MJD 
54682 - 55414). We used the standard LAT analysis soft- 
ware, ScienceTools v9r21. The events were selected using 
so-called "diffuse class" events. In addition, we excluded 
the events with zenith angles greater than 100° to avoid 
the contamination of the Earth-limb secondary 7 radi- 
ation. The events were extracted in the range between 
200 MeV and 300 GeV within a 15° acceptance cone of 
the Region of Interest (ROI) centered on the location of 
3C 279 (RA = 195.047°, DEC=-5.789°, J2000). Below 
200 MeV, the effective collection area of LAT for the dif- 
fuse class events drops very quickly and thus larger sys- 
tematic errors are expected. The 7-ray flux and spectrum 
were calculated using the instrument response function 
(IRF) of "P6_V11_DIFFUSE" by an unbinned maximum 
likelihood fit of model parameters. We examined the sig- 
nificance of the 7-ray signal from the sources by means 
of the test statistic (TS) based on the likelihood ratio 
tes10 The background models included a component 
for the Galactic diffuse emission along the plane of the 
Milky Way, which was modeled by the map cube file 
"glLierrk.v02-P6-Vll-DIFFUSE.fiti? . An isotropic com- 
ponent (isotropic-iem-v02-P6-Vll-DIFFUSE.txt) was 

59 TS = 25 with 2 degrees of freedom corresponds to an es- 
timated ~ 4.6 a pre-trials statistical significance assuming that 
the null-hypothesis TS distribution follows a x 2 distribution (see 
IMattox et aT1[T996Tl . 
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also included to represent the extragalactic diffuse emis- 
sion and residual instrumental background. Besides 
those components, the model in our analysis also in- 
cluded the emission from all nearby point sources in- 
side the ROI from the first Fermi-LAT catalog (1FGL: 
lAbdo et al.ll2010bt ). The spectra of those sources were 
modeled by power law functions except for a pulsar 
1FGL J1231. 1-1410 (=PSR J1231-1411), for which we 
included an additional exponential cut-off in its spectral 
modeling. During the spectral fitting, the normalization 
factors of the Galactic diffuse and isotropic components 
and the nearby sources were left as free parameters, and 
the photon indices of the nearby sources were fixed to the 
values from the 1FGL catalog except for 3C 273, whose 
photon index was allowed to vary freely. In the light 
curve analysis, we considered only two bright sources in 
the background model as nearby point sources, namely 
3C 273 and 1FGL J1231. 1-1410, because other nearby 
sources had a negligible contribution to 7-ray signal, es- 
pecially in such relatively short time scales (shorter than 
a week) for the light curves considered here. The fluxes 
used for the light curve were calculated by a simple power 
law model fit using data in the given energy ranges. 



2.2. Temporal behavior 

The 7-ray light curve measured by Fermi-LAT can be 
seen in Figure[T] The Figure shows the flux history above 
200 MeV averaged over (a) 1-day intervals, (b) 3-day in- 
tervals, and (c) 1-week intervals. It also includes 1-week 
light curves of (d) the flux between 200 MeV and 1 GeV, 
(e) the flux above 1 GeV, and (f) the photon index in 
the range above 200 MeV. 

The 7-ray flux clearly shows variability. The source 
showed high-flux states between MJD 54700 and 54900, 
in which two prominent flares can be seen: one of 
the flares at ~ MJD 54800 and the the other at ~ 
MJD 54880. During the second flare, a change in the 
optical polarization associated with a 7-ray flare was 
discovered (Paper I). We detected some flux variabil- 
ity between MJD 55000 and 55120, but after that, the 
source remained in a relatively low-activity state until 
the end of the period considered in this paper. Dur- 
ing this 2-year period, the highest integral flux above 
200 MeV occured on MJD 54880 in the 1-day interval 
light curve with flux of Fe >2 oomcV= (11-8 ± 1.5) x 10~ 7 
photons cm~ 2 s _1 and TS — 306. By extrapolating 
the spectrum down to 100 MeV, an integral flux above 
100 MeV on that day yields F B >iooMeV= (31.0 ± 6.0) x 
10~ 7 photons cm~ 2 s _1 , which is still a factor of 3-4 
times lower than the flux of the brightest flare (~ 1 x 10~ 5 
photons cm -2 s _1 ) d etected during the EGRET obser- 
vation s of the source (|Wehrle et al.lll998l : lHartman et al.l 
I2001bl) . 

We quantified the flux variability using 1-week inter- 
val data for energies above 200 MeV [full band] , between 
200 MeV and 1 GeV [soft band], and above 1 GeV [hard 
band]. This is based on the "excess va riance" method 
(|Nandra et al.lll997t lEdelson et al.ll2002ft after subtract- 
ing the contribution expected from measurement errors 
(ccnyi)- Using the mean square error < < 7 e rr.i >. the ex- 
cess variance F var can be described as (jVaughan et al.l 



12001 

_ I S*-<* c „,i> [ 

P ™ - V < F > 2 1 ' 

where S is the variance of the flux, and < F > is the 
mean value of the flux. The definition o f associated error 
can be found in (jVaughan et al.l |2003|) . In the calcula- 
tion, we excluded bins of 8, 82, 87 and 90 because the 
fit in the hard band failed due to poor statistics of the 
data samples. Resulting F var values are 0.695 ± 0.015, 
0.648 ±0.017 and 0.839 ±0.030 for the full, soft and hard 
bands, respectively. The resulting values indicate that 
the flux of the hard band showed significantly stronger 
variability than that of the soft band. For comparison, 
F var = 0.79 ± 0.02 for E > 300 MeV has been reported 
durin g the first 11 mon ths of the Fermi scientific mis- 
sion ([Abdo et al.1 l2010fl ) , when the source has clearly 
been more active. 

A Power Density Spectrum (PDS) for the 3-day binned 
light curve was calculated using a Fourier transform 
and is shown in Figure [2j The power density was 
normalized to fractional variance per frequency unit 
( rms 2 I~ 2 Day^ 1 ) and the PDS points were averaged in 
logarithmic frequency bins. The white noise level was 
estimated from the rms of the flux errors and was sub- 
tracted from the PDS. A slope of 1.6 ± 0.2 was obtained 
from a linear fit to the binned PDS for frequencies up 
to 0.1 Day -1 . The main uncertainty in the estimated 
PDS slope is due to the stochastic nature of the variabil- 
ity which leads to variations in the determined slope be- 
tween different time limited observations. An additional 
effect which can cause a systematic bia s in the observed 
PDS slope is the red noise leakage (e.g. JChatteriee et al.l 
2008) . In the present analysis this effect is not taken into 
account. 

Figure [3] shows plots of flux vs. photon index (F) 
based on the weekly light curve results above 200 MeV 
[full band], between 200 MeV and 1 GeV [soft band] 
and above 1 GeV [hard band] . The data which have 
TS > 10 were selected for the plots and are shown in 
gray points. An average photon index was calculated by 
fitting a constant value in each plot, corresponding to 
r>200MeV = 2.334 ± 0.015, r 2 ooMoV-iGcV = 2.20 ± 0.03 
and r > i GoV = 2.48 ± 0.04 for the full, soft and hard 
bands, respectively. The average photon index in the soft 
band shows a significantly harder spectrum than that in 
the hard band. 

We also derived photon indices resulting from an 
analysis where the data were sorted in five bins using 
week-long fluxes for each energy band, and plotted the 
results as red points. Those photon indices of each flux 
bin are also shown in the insets in Figure |3] For the 
full band, although the change of the photon index is 
rather small (Ar ~ 0.2) compared to the flux variation 
(spanning about an order of magnitude) , a weak "harder 
when brighter" effect can be seen. Suc h an effect 
was al so measured in other LAT blazars ()Abdo et al.l 
l2010dt l. The soft band also shows the weak "harder 
when brighter" effect with a slightly larger change of 
the photon index (Ar ~ 0.4). On the other hand, the 
photon index of the hard band changes only sightly 
(Ar ~ 0.1) and is statistically consistent with a constant 
value. 
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Figure 1. Gamma-ray light curves of 3C 279 during the first two years of the Fermi-LAT observations from 2008 August to 2010 August, 
plotted in; (a) 1-day intervals at energies above 200 MeV, (b) 3-day intervals at energies above 200 MeV, (c) 1-week intervals at energies 
above 200 MeV, (d) 1-week intervals at energies between 200 MeV and 1 GeV, (e) 1-week intervals at energies above 1 GeV. The panel- (f) 
shows the history of the photon index at energies above 200 MeV in 1 week intervals, while the panel- (g) shows arrival time distribution 
of > 20 GeV events associated with 3C 279. The vertical axis of the panel-(g) represents the estimated energy of events. The highest 
energy photon corresponds to 30.8 GeV at MJD 54891. The dotted lines and capital letters represent time intervals where 7-ray spectra 
are extracted (see also Table [TJ- 

2.3. Highest energy photons 

During the 2-year observations, the highest en- 
ergy photon associated with 3C 279 was detected at 
MJD 54891.60745 with an estimated energy of 30.8 GeV. 
The event was converted in the front-thin layers (so- 
called "front event") of the LAT detector and still re- 
mains even when we apply the cleanest event selec- 
tion (so-called "data clean event" ) , which w as developed 
for stu dying extragalactic 7-ray background (|Abdo et al.l 
2010g). The reconstructed arrival direction of the event 
is 5'. 7 (=0°.095) away from 3C 279, and is within the 68% 
containment radius of the LAT PSF (0°.114 in the IRF 
of "P6.V11.DIFFUSE") for front events at 30.8 GeV. 
Based on our model fit of the epoch which contains that 
highest-energy photon, we find the probability that the 
photon was associated with 3C 279 (as opposed to all 
other sources in the model including the diffuse emission 
and nearby point sources) is 88.6 %. 

In total, we found 10 events with estimated energies 
higher than 20 GeV within an 0°.25 radius centered at 
3C 279. All events lay within a 95% containment radius 



CD 
O 




-2.0 -1.5 
LOG FREQUENCY (day 1 ) 



0.5 



Figure 2. Power Density Spectrum of 3C 279 for the 3-day binned 
7-ray light curve. The white noise level has been subtracted. The 
solid line histogram describes the PDS averaged in logarithmic fre- 
quency bins while the dotted curve describes the raw PDS before 
binning. The dashed line represents a linear fit to the binned PDS. 
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Figure 3. Plots of the integrated 7-ray flux vs. photon index of 3C279 measured in week-long bins for energies above 200 MeV (A), 
between 200 MeV and 1 GeV (B) and above 1 GeV (C). Only points with TS > 10 are plotted (gray points). The blue dotted horizontal 
lines indicate average photon indices of those data for each energy band. Red points show the photon indices resulting from an analysis 



where the data were sorted in five bins using week-long fluxes for each energy band, 
ranges of the week-long flux bins while the vertical bars indicate 1 a statistical errors, 
well as the average photon index of each energy band. 



For the red points, the horizontal bars indicate the 
The insets show enlarged views of the red points as 



of the LAT PSF from 3C 279 and remain even after the 
"data clean selection" applied. The number of expected 
background events above 20 GeV within the 0°.25 radius 
at the location of 3C 279 for the 2-year observations is 
only 0.16 events. The bottom panel of Figure Q] plots the 
arrival time distribution of those 10 events. All events 
except for two were detected between MJD 54780 and 
54900 during the high activity states. No photon above 
20 GeV associated with 3C 279 has been detected after 
MJD 54914 during the 2-year observations. 

2.4. Gamma-ray spectra 

We extracted the 7-ray spectra using data for the en- 
tire 2-year period and following 8 sub-periods (see also 
Table [l}: [A] the initial quiescent state in the 7-ray 
band (MJD 54682 - 54728), [B] the first 7-ray flaring 
state (MJD 54789 - 54809), [C] an intermediate state 
(MJD 54827 - 54877), [D] the first 5 days of the sec- 
ond 7-ray flaring event (MJD 54880 -54885), [E] the last 
3 days of the second 7-ray flaring event (MJD 54897 - 
54900), [F] during the isolated (first) X-ray flaring event 
(MJD 54950 - 54960; see Section HD, [G] during the sec- 
ond X-ray flaring event (MJD 55042 - 55045; see Sec- 
tion EJ), and [H] a quiescent state (MJD 55240 - 55319). 
Those sub-periods were also selected taking into account 
observations in other energy bands. Spectral energy dis- 
tributions in the 7-ray band for each sub-period are pre- 
sented in Figure [U Each 7-ray spectrum was modeled 
using a simple power law (PL; dN/dE oc E~ r ), a bro- 
ken power law (BPL; dN/dE oc E~ Tl for E < E hl± and 
dN/dE oc E~ T2 otherwise), and a log parabola (LogP; 

dN/dE oc (E/E )~ a ~ p] ° e( - E/Eo) ) model. In the case of 
LogP model, the parameter /3 represents the curvature 
around the peak. We note that the choice of the refer- 
ence energy E$ in the LogP model does not affect the 
determination of the other two model parameters, and 
hence we fixed it at 300 MeV. 

The best-fit parameters calculated by the fitting pro- 
cedure are summarized in Table [1] The integral fluxes 
above lOOMe'VFl derived using each spectral model are 
also included. The averaged 7-ray spectral shape for 

60 Although we use photon data from 200 MeV, the integral 
fluxes are extrapolated down to 100 MeV, which is convenient to 
compare with other 7-ray results. 
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Figure 4. Gamma-ray spectral energy distributions of 3C 279 
of each period as defined in the text or Table [l] 2-year averaged 
[black: filled circles], Period A (magenta: open circles), Period B 
(brown:f illed triangles), Period C (orange: open triangles), Period 
D (red: filled squares), Period E (green: open squares), Period F 
(pink: filled diamonds), Period G (cyan: open diamonds) and Period 
H (blue: crosses). The vertical bars indicate 1 <r statistical errors 
while the horizontal bars indicate energy ranges of each bin. The 
best-fit spectral models are plotted as dotted lines for each period 
and their parameters are summarized in Table[l] We use the broken 
power law model for the spectra of 2-year, Period B and Period C 
because significant improvements in the spectral fits can be seen 
compared to the simple power law model (see Table [TJ while the 
simple power law model is used for other periods. The lower panel 
shows the residuals, plotted as X (= (data — model)/data error) 
from the best-fit models. A 'dip' feature at ~ 1—2 GeV in the 
spectrum of Period G (the 3rd point in cyan) is a < 2 <r effect from 
the best-fit model, thus, not statistically significant. 
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Table 1 

Results of spectral fitting in the 7-ray band measured by Fermi-LAT. 



Period 
(MJD) 


fitting niodcl a 


Gamma-ray spectrum (Fermi-LAT) 
r/a/ri /3/r 2 E hlk (GeV) 


TS 


-2AL b 


Flux (> 100 MeV) 
(10 -7 ph cm -2 s _1 ) 




PL 


2.38 ±0.02 






20272 




U . 1 U _1_ U.1J 


ZUUo -rt.Ug 4 — ZU1U -A-Ug 


T nrrT3 
LOgf 


2.18 ± 0.03 


0.08 ± 0.01 




ZUZO t 


AP. c; 
40. 


O.lo it U.10 


|^D40oZ — D0414J 


DDT 

DE L 


2.31 ±0.02 


2.95 ±0.12 


3.5 ±0.3 


ZUZoD 


A1 O 
4o.U 


x. TO, 4-f| K 
O. (D ± U.lO 


fcriod A 


PL 


2.30 ±0.07 






797 




3.7 ± 0.4 


ZUUO -rVUg 4 — ZUUo oep iy 




2.19 ± 0.15 


0.04 ± 0.05 






n 1 


o.o ± U.O 


(54682 - 54728) 


BPL 


2.21 ±0.09 


2.82 ± 0.40 


3.4 ±0.8 


798 


2.5 


3.4 ±0.4 


± Cl 1UU _L> 


PL 


2.28 ±0.04 






3209 




19.0 ± 1.1 


9008 Mnv 1 Q 900S T\ar* Q 

zuuo r\ov iy — zuuo uec y 


T ncrP 


1.95 ± 0.10 


0.13 ± 0.04 




391 A 


lO.O 


±o.u it i.o 


(54789 - 54809) 


BPL 


2.00 ±0.10 


2.61 ±0.11 


1.0 ±0.2 


3215 


13.7 


15.8 ± 1.3 


± <_ 1 ll_Jt_l. 


PL 


2.25 ±0.04 






4107 




i n n + o 5 

1U.U _1_ W. iJ 


zuuo ucc z; — zuuy reo 10 




2.05 ± 0.08 


0.08 ± 0.03 




4±1U 


o. Z 


o.O it U.D 


(54827 - 54877) 


BPL 


2.07 ±0.08 


2.43 ±0.08 


1.0 ±0.2 


4109 


8.2 


8.9 ±0.6 


Period 13 


PL 


2.36 ±0.08 






1236 




23 6 + 2 3 


zuuy pcd io — zuuy reD zo 




2 1 6 ± 1 6 


n no, + D6 




1 OQ/1 
1Zo4 


z.u 


on 9 _l_ 9 c 
ZU.Z it Z.o 


^04ooU — 04oo0^ 


DDT 


2.25 ±0.12 


2.91 ±0.61 


2.3 ±2.1 


IZoO 


Z.y 


Ol Q 1 O P 

Zi.o ± Z.O 


Period E 


PL 


2.64 ±0.32 






61 




6.3 ±2.5 


2009 Mar 7 - 2009 Mar 10 
(54897 - 54900) 


LogP 


2.64 ±0.32 


0.00 ±0.00 




61 


0.0 


6.3 ±2.5 


Period E 


PL 


2.54 ±0.24 






85 




3.5 ± 1.2 


2009 Apr 29 - 2009 May 9 
(54950 - 54960) 


LogP 


2.54 ±0.24 


0.00 ±0.00 




85 


0.0 


3.5 ± 1.2 


Period G 


PL 


2.44 ±0.13 






460 




18.8 ±2.9 


2009 Jul 30 - 2009 Aug 2 
(55042 - 55045) 


LogP 


2.37 ±0.25 


0.03 ±0.10 




460 


0.1 


17.7 ±4.0 


Period H 


PL 


2.83 ±0.11 






398 




3.7 ±0.5 


2010 Feb 13 - 2010 May 3 


LogP 


2.56 ±0.23 


0.16 ±0.13 




399 


1.9 


2.9 ±0.6 


(55240 - 55319) 


BPL 


2.72 ± 0.43 


3.47 ±0.57 


1.6 ±0.5 


399 


1.6 


3.4 ± 1.7 



PL: 

AL represents the diffcrc 



power law model, LogP: log parabola model, BPL: broken power law model. Sec definitions in the text. 

: of the logarithm of the likelihood of the fit with respect to a single power law fit. 



the 2-year observation significantly deviates from a sin- 
gle power law. A LogP model is favored to describe 
the 7-ray spectral shape over the simple PL model with 
the difference of the logarithm of the likelihood fitj^l 
—2AL = 46.5 (corresponding to a significance level of 
- 7(70 and a BPL fit yields -2AL = 43.0. Even 
in some individual periods as defined above, the spec- 
tra deviate from a single power law: for example, the 
spectrum in the Period B yields —2AL = 13.6. This 
is consistent with our finding in Section 2.2 that the 
spectrum above 1 GeV is significantly softer than the 
spectrum below 1 GeV. We thus conclude that the 7-ray 
spectrum significantly deviates from a simple power law. 
The spectral break in 3C 279 is not as p ronounced as that 
seen in the spectra e.g., of 3C 454.3 ()Ackermann et al.l 
2010). One the other hand, the BPL model returns 
break energies within a few GeV range regardless of the 
flux levels as observed in other bright FSRQ s, such as 
3C 454.3 and 4C+21.35 (jTanaka et al.l 1201 if ). Such a 
spectral feature could be due to 7 — 7 absorption to 
pair production by He II Lyman recombination contin- 
uum UV photons from t he emission line region (see, e.g., 
iPoutanen fc Sternl[20lol). or a break in the electron dis- 
tribution (|Abdo et alJl2009bD . We consider the 7-ray 
emission region to be located significantly beyond the 
broad emission region (see the Discussion in Section [5]), 
and this implies that the break in the electron energy 
distribution is the more likely explanation. 

61 -2AL = -21og(L0/Ll), where L0 and LI are the maxi- 
mum likelihood estimated for the null and alternative hypothesis, 
respectively. 

62 Because the LogP model has one more free parameter than 
the PL model has, the — 2AL distribution follows a \ 2 distribution 
with 1 degree of freedom. 



3. MULTI- WAVELENGTH OBSERVATIONS AND DATA 
REDUCTION 

3.1. X-ray and Hard X-ray: Suzaku 

The Suzaku X-ray satellite (|Mitsuda et al.l I2007T ) ob- 
served 3C 279 as a part of multi-band studies of the 
object. The observations took place in two segments, 
with an interruption lasting roughly 1.5 days: (1) be- 
tween 2009 January 19, 23:19:00 and 2009 January 22, 
22:32:00 UTC (sequence number 703049010), and (2) be- 
tween 2009 January 23, 20:45:00 and 2009 January 25, 
03:00:00 UTC (sequence number 703049020). "Period 
C" (see Tabled]) includes both Suzaku observations. The 
goals of the Suzaku observations were to monitor the soft- 
medium X-ray flux (0.3-12 keV) o f the source with the 
X-ray Imaging Spectrometer (XIS: iKoyama et all 120071 ) 
and to take advant age of the data from the Hard X- 
ray Detector (HXD; iTakahashi et al.ll2007t ). The HXD 
consists of PIN silicon diodes for the lower energy band 
(10-70 keV) and GSO scintillators for the higher energy 
band (40-600 keV), to extend the spectral bandpass be- 
yond the energies accessible with imaging instruments 
(> lOkeV). The HXD nominal position was used for the 
observations to maximize its effective area. In the follow- 
ing analysis, the HXD/GSO data were not used because 
there was no significant detection of the source. 

Although the observation conditions were nominal, the 
XIS1 data suffered from somewhat high and variable 
background, resulting in the total apparent counting rate 
ranging from 1 to 3 counts s _1 in source- free regions for 
the entire chip. Still, the background-subtracted spec- 
trum determined from the XIS1 data below 8 keV was 
entirely consistent with that from XIS0 and XIS3 and 
thus we included the background-subtracted XIS1 data 
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in the spectral fitting. The total duration of good data 
accumulated by the XIS instruments was 191 ks. We 
used the standard f tools data reduction package, pro- 
vided by the Suzaku Science Operations Center, with 
the calibration files included in CALDB ver. 4.3.1. For 
the analysis of spectra and light curves, we extracted 
the counts from a region corresponding to a circle with 
260"radius, centered on the X-ray centroid; we used a 
region of a comparable size from the same chip to ex- 
tract the background counts. The net count rates were 
0.47, 0.63, and 0.56 count s" 1 for XISO, XIS1, and XIS3, 
respectively, with the typical count rate uncertainty in 
the entire observations of ~ 3 %. The data indicate no 
significant variability during the Suzaku observations. 

The source was also detected in the HXD/PIN data, al- 
though the signal was relatively weak. We used the stan- 
dard cleaned events, processed using the standard crite- 
ria applicable to the rev. 2.13 of the Suzaku HXD data 
processing software. This yielded 95.4 ks of good data, 
with a total count rate of 0.3 count s _1 . For the back- 
ground subtraction, we used the standard background 
files provided by the Suzaku team through HE AS ARC. 
We applied the standard tool hxdpinxbpi which ac- 
counts for the particle background as well as for the con- 
tribution of the Cosmic X-ray Background as appropri- 
ate for the effective area and the solid angle of the HXD. 
The net counting rate was 0.02 count s _1 , with the for- 
mal statistical uncertainty of ~ 10%. We note that this 
formal uncertainty is probably lower than the standard 
systematic error due to the background subtraction of 
3% of the average background (corresponding to 0.01 
count s _1 ). Nonetheless, even if the additional uncer- 
tainty is included, the source was still detected by the 
HXD/PIN. 

For the spectral analysis, we used the XSPEC spectral 
analysis software. For the spectral fitting of the XIS 
data, we used the standard redistribution files and mir- 
ror effective areas generated with Suzaku-specific tools 
xisrmf gen and xissimarf gen. In the spectral fits, we 
used the counts corresponding to the energy range of 
0.5-10.0 keV for XISO and XIS3, and 0.5-8.0 keV for the 
XIS1. We used all three XIS detectors simultaneously, 
but allowed for a small (a few %) variation of normaliza- 
tion. For the HXD/PIN data, we considered the data 
in the range of 20-50 keV and used the response file 
ae_tixd_pinhxnome5_20080716 . rsp. 

The source spectrum was modeled as an absorbed 
power law, with t he cross-sections and elementa l abun- 
dances as given in lMorrison fc McCammonl (|1983[ ): other 
absorption models give similar results. The best-fit ab- 
sorbing column was (3.1 ± 0.5) x 10 20 cm~ 2 , and the 
photon index was 1.76±0.01. Inclusion of the HXD/PIN 
data in the fit did not change the fit parameters percep- 
tibly. The x 2 f° r the fit including the three XIS de- 
tectors and the HXD/PIN was acceptable, with 5061 
for 5023 PHA bins. The absorption inferred from the 
simple absorbed power law model is marginally greater 
than the value inferred from the radio measurements of 
the column density of the material in the Galaxy of 
2.0 x 10 20 cm" 2 (wi th an estimated error of ~ 10%) 
(|Kalberla et al.ll2005l ). We deem the difference not sig- 
nificant, since at such small column densities, it can be 
accounted for by even small systematic uncertainty in 



the knowledge of the effective area of the XIS instru- 
ments at the lowest end of the XIS bandpass. Further- 
more, a modest additional column density is expected 
in the host galaxy of 3C279. The observed model 2-10 
keV flux is 8.0 x 10~ 12 erg cm~ 2 s _1 , with the statistical 
error of < 2%, which is probably smaller than the sys- 
tematic error resulting form the calibration uncertainty 
of the Suzaku instruments. We plot the Suzaku 3C 279 
spectra in Figure [5j 
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Figure 5. Count spectra of 3C 279 measured by Suzaku XISO 
(black), XIS 1 (red), XIS3 (green) and HXD/PIN (blue). The model 
plotted with the data is a broken power law obtained by a fitting 
these three XISs and HXD/PIN data. The lower panel shows the 
residuals for this broken power law model. 

3.2. X-ray: XMM-Newton 

XMM-Newton observed 3C 279 once starting on 2009 
January 21, 17:28 UT. The observation was largely de- 
void of flares (except for the period close to the end of the 
observation), and the total length of good data accumu- 
lated in the pointing was 16.8 ks. We used the standard 
Scientific Analysis System (SAS) data reduction pack- 
age, provided by the XMM-Newton Science Operations 
Center. Since 3C 279 is a relatively bright source, we 
considered only the pn-CCD data. We note here that 
the spectra and light curves taken by MOS-CCDs were 
entirely consistent with the results inferred from the pn- 
CCD data. 

For the analysis of spectra and light curves, we ex- 
tracted the counts from within 40"radius of the source; 
we used a region of the same size, from the same pn-CCD 
chip, to extract the background counts. The data indi- 
cate no significant variability during the XMM-Newton 
observation. The spectral analysis was performed us- 
ing the XSPEC v.12 spectral analysis software with the 
standard redistribution files and mirror effective areas 
included in the SAS package. We used the counts cor- 
responding to the energy range of 0.5-10.0 keV in our 
spectral fits. 

The source spectrum was first modeled as an ab- 
sorbed power law; the best-fit absorbing column was 
(2.2 ± 0.6) x 10 20 cm" 2 , and the photon index was 
1.77 ± 0.03, with x 2 of 588 for 517 d.o.f. The result is 
consistent with the spectral analysis results of the Suzaku 
observations as described in the previous section, which 
were performed during the same period as the XMM- 
Newton observation. We also considered a broken power 
law model and found that the overall intrinsic source 
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spectrum hardens with increasing energy. The absorb- 
ing column was (3.4 ± 0.7) x 10 20 cm -2 , and the low 
and high energy indices were respectively 1.83 ±0.05 and 
1.55 ± 0.2 with the break energy of 4.1 ± 0.8 keV. The 
resulting y 2 was 563, for 515 d.o.f. The broken power 
law model is statistically only marginally superior to the 
simple power law model, especially given that the absorp- 
tion inferred form the s imple power law mode l is closer to 
the value inferred from lKalberla et al.l (|2005f) . For either 
model, the 2-10keV flux is 7.7 x 10" 11 erg cm" 2 s" 1 , 
with a statistical error of 5%, which is probably smaller 
than the systematic error resulting form the calibration 
uncertainty of the XMM- Newton pn-CCD. 

3.3. X-ray: RXTE-PCA 

RXTE carried out 321 observations between 2008 July 
3rd (MJD 54650) and 2010 August 12th (MJD 55420). 
Those include 52 observations based on the Cycle 12 
Guest observer (GO) program and 269 observations 
based on the Core program in Cycle 12-14. The fluxes 
resulting from the Cyclel2 GO observations have been 
already reported in Paper I. Most of the individual obser- 
vations have exposure times in a range from 1.0 to 2.5 ks. 
We analyzed the data from the Proportional Counter 
Array (PCA) following standard procedures using the 
rex script in HEASOFT v. 6. 9. Only signals from the 
top layer (X1L and X1R) of PCU2 were extracted for 
data analysis. The data were screened with the follow- 
ing data selection: source elevation above the horizon 
> 10°, pointing offset smaller than 0.02°, at least 30 
minutes away from a SAA passage and electron contam- 
ination smaller than 0.1. The background was estimated 
with standard procedures, and the detector response ma- 
trices were extracted with the RXTE tools (command 
PCARSO v. 11. 7.1). For the spectral analysis we re- 
binned the spectra into 11 channels. The spectra from 
the channels corresponding to nominal energies of 2.6 
to 10.5 keV are adequately fitted by a single power law 
model, absorbed by a fixed Galactic column density of 
2.2 x 10 20 cm -2 using the XSPEC v. 12 software package. 
The value of the column density is based on our XMM- 
Newton results (in Sectio nEL2J and is also consistent with 
the value based on iKalberla et al.l (|2005l ). 

3.4. X-ray: Swift-XRT 

In the HEASARC data baseQ there are 80 pub- 
licly available Swift X-Ray Telescope (XRT) observa- 
tions between 2008 July 3 (MJD 54650) and 2010 Au- 
gust 12 (MJD 55420), which include 32 pointings based 
on an approved GI proposal in Cycle-4 (Proposal num- 
ber: 5080069). The results of the flux history based 
on the data until 2009 May 31 have already reported 
in Paper I. Effective exposure times of these observa- 
tions range between 1 and 3 ks, but some have longer 
exposure times, for example, 8.9 ks for ID:35019007 
(MJD 54795), 22.5 ks for ID:35019009 (MJD 54797), 
20.2 ks for ID:35019010 (MJD 54799) and 15.4 ks for 
ID:35019011 (MJD 54800). The XRT was used in the 
photon counting mode, and no evidence of pile-up was 
found. The XRT data were reduced with the standard 
software xrtpipeline v. 0.12.6, applying the default fil- 
tering and screening criteria (HEADAS package, v. 6. 10). 

6,5 http : //heasarc .gsf c .nasa. gov/cgi-bin/W3Browse/ swift .pi 



The source events were extracted from a circular region, 
20 pixels in radius, centered on the source position. Ex- 
posure maps were used to account for Point-Spread Func- 
tion losses and the presence of dead pixels/columns. The 
background was determined using data extracted from 
a circular region, 40 pixels in radius, centered on (RA, 
Dec: J2000) =(12h56m26s,-05°49'30"), where no X-ray 
sources are found. Note that the background contamina- 
tion is less than 1 % of source flux even in the faint X-ray 
states of the source. The data were rebinned to have at 
least 25 counts per bin, and the spectral fitting was per- 
formed using the energy range between 0.3 keV and 10 
keV using XSPEC v. 12. The Galactic column density is 
fixed at 2.2 x 10 20 cm -2 during the fittings as is the case 
in the RXTE data analysis. 

Figure[6]shows a scatter plot between photon index and 
flux in the X-ray band as measured by Swift-XRT and 
RXTE-PCA. Generally, a "harder-when-brighter" trend 
can be seen. Only the highest-flux point measured by 
Swift-XRT shows the photon index significantly harder 
(smaller) than 1.5. 
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X-ray flux (2-10keV) 



Figure 6. Scatter plot of flux vs. photon index of 3C 279 in the X- 
ray band with the data taken by Swift-XRT (magenta) and RXTE- 
PCA (cyan). The horizontal dotted line represents the photon 
index value of 1.5. 

3.5. Ultra-Violet: Swift-UVOT 

The Swift Ultra-Violet/Optical Telescope (UVOT; 
Roming et al. 2005) data used in this analysis included 
all of the observations performed during the time in- 
terval MJD 54650-55420. The UVOT telescope cycled 
through each of the six optical and ultraviolet filters 
(V,B,U,W1,M2, W2). The UVOT photometric sys- 
tem is described in lPoole et al.l (|2008l ). Photometry was 
computed from a 5" source region around 3C 279 us- 
ing the publicly available UVOT FT00LS data reduction 
suite. The background region was taken from an annu- 
lus with inner and outer radii of 27". 5 and 35", respec- 
tively. Galactic absorpt ion in the dire c tion o f 3C 279 
was adapted as given in lLarionov et al.l (|2008f ) , namely, 
A v = 0.093, A B = 0.123, A v = 0.147, A W1 = 0.195, 
Am2 — 0.285 and Aw2 — 0.271. The measured mag- 
nitudes in each band during the 2-year observations are 
m v = 15.6 - 18.7 (75 data points), m B = 16.0 - 18.8 
(80 data points), mjj = 15.1 — 18.0 (88 data points), 
mwi = 15.3 - 18.0 (84 data points), m M 2 = 15.4 - 18.2 
(76 data points) and mw2 = 15.5—18.0 (81 data points). 
All observed data points are shown in Figure [Jj 
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Figure 7. Light curves of all observed UV-optical-NcarlR bands 
of 3C 279 in measured magnitude scale from 2008 August to 
2010 August, including W2 (Swift-VVOT), M2 (Swift-VVOT), Wl 
(Swift-VVOT), U (Swift-VVOT), B (Swift-VVOT), V (Katana, 
Swift-VVOT), R (Abastumani, Calar Alto, ST-7, GRT, MDM, 
L'Ampolla, Perkins, SLT, KVA, LT, San Pedro, St. Petersburg, 
Tijarafc), J (AZT-24, Kanata), H (AZT-24, LT) and K (AZT-24, 
Kanata) bands. 

3.6. Optical, Near- Infrared and Radio observations by 
GASP-WEBT 

The GLAST-AGILE Support Program (GASP; 
IVillata et all 120081 120091 ) is a project initially originat- 
ing from the Whole Earth Blazar Telescop<H (WEBT) 
in 2007. It is aimed to provide long-term monitoring in 
the optical (i? band), near-IR, and mm-cm radio bands 
of 28 7-ray-loud blazars during the lifetime of the AGILE 
and Fermi 7-ray satellites. 

The observations of 3C 279 in the period considered 
in this paper were performed by the observatories listed 
in Table [2] The calibrated i?-band magnitudes of the 
source were obtained through differential photometry 
with respect to th e reference stars 1, 2, 3, and 5 by 
iRaiteri et ail (|1998[ ). Near-IR data in the J, H, and K 
filters were acquired at Campo Imperatore and Roque 
de los Muchachos (Liverpool). When converting magni- 
tudes into flux densities, optical and near-IR data were 
corrected for Galacti c reddening using Ab = 0.123 mag 
(Schl egeTet al.l 119981). W e adapted the extinction laws 
by Cardelli et al.l~j 1989f) . and the zero-mag fluxes by 
IBessell et all (11998ft . 

For the observations between 2008 August and 2010 
August, the measured i?-band magnitude ranged from 
14.87 to 17.81 (673 data points). The i?-band data have 
the best time coverage among the IR-optical-UV bands in 
our data thanks to the participation of a number of tele- 
scopes. The emission shows strong variability and the ex- 

64 http : //www. oato . inaf . it/blazars/webt 



Figure 8. Radio light curves of 3C 279 from 2008 August to 
2010 August measured at all observed radio frequencies: 345 GHz 
(SMA), 230 GHz (CARMA, SMA), 93 GHz (CARMA), 43 GHz 
(Noto), 37 GHz (Metsahovi), 22 GHz (Medicina), 15 GHz (OVRO), 
14.5 GHz (UMRAO), 8 GHz (Medicina, UMRAO) and 5 GHz 
(Medicina, UMRAO). 

cess variance (F v&1 ; Eq.[l| of the source i?-band flux (i.e., 
in linear scale) is 0.853 ± 0.001. The near-IR magnitudes 
in the J, H and K bands were measured in ranges of 
m, } = 14.91 - 15.59 (20 data points), m H = 12.04-14.90 
(68 data points) and mx = 11.19 — 13.45 (20 data 
points) . Those data points are shown in Figure [7J The 
radio flux densities were measured in ranges of F5 ghz = 
8.5-12.4 Jy (109 data points), F 8GHz = 9.1-15.5 Jy (124 
data points), Fu 5 ghz = 10.3 — 19.4 Jy (118 data points), 
F22 GHz = 10-22 Jy (16 data points), F 37GH z = 10-20 
Jy (168 data points), F 43 gh z = 10-22 Jy (20 data 
points), i*230 GHz = 5.1 — 10.5 Jy (62 data points) and 
-F345GH Z = 6.0 — 6.8 Jy (7 data points). The light curves 
of the radio flux densities in those bands are plotted in 
Figure M 

3.7. Optical and Near- Infrared: the Kanata telescope 

We performed the V, J and Xs-band photometry and 
polarimetry of 3C 279 using TRISPEC installed to the 
1.5 m Kanata telescope located in the Higashi- Hiroshima 
Observatory. 

TRISPEC has a CCD and two InSb arrays, enabling 
photopolarimetric observations in an optical and tw o 
near-IR bands simultaneously ( Watana be et al.l |2005|) . 
We obtained 64, 42 and 17 photometric measurements 
in the V, J and Ks bands, respectively. A unit of 
the polarimetric observing sequence consisted of suc- 
cessive exposures at 4 position angles of a half-wave 
plates: 0°, 45°, 22°. 5, 67°.5. The data were reduced ac- 
cording to the standard procedures of CCD photometry. 
We measured the magnitudes of objects with the aper- 
ture photometry technique. We performed differential 
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Table 2 

List of observatories that contributed data to the campaign 



O bscrvatory 


Detector/ telescope (diam.) 


Band 




Oamrna ray 




1 O / 1 ILL 


I j . 1 1 1 BUI V \_- y ±11UCH_. 1 


■> 200 MpV 




X ray 




Suzaku 


XIS 


0.5 -10 keV 




HXD/PIN 


15 - 50 keV 


XMM-Newton 


PN 


0.5 -10 keV 


RXTE 


PCA 


3-10 keV 


Swift 


XRT 


0.6 - 7 keV 




Ultra- Violet, Optical, Infrared 




Swift 


UVOT 


W2,M2, Wl, U, B, V 


Spitzer 


IRS 


5-38 /an 




TT) (VP 


O.O, ^t.O, O.o, O.U fllll 




MIPS 


24, 70, 160 /im 


Abastumani, Georgia* 


(70 cm) 


R 


Calar Alto*'* 






Campo Imperatore, Italy 1 ^ 


AZT-24 (110 cm) 


J, H, K 


Crimean, Ukraine^ 


ST-7 (70 cm) 


R 


p A j,i „,,,] TTQA# 


\jsx\- 


IX 


Hiroshima, Japan 


Kanata (150 cm) 


V, J, Ks, polarization (V) 


1Y1 L L JT t_ ctiv , ill IZjUllel, U On 


MT)A/T (1 3D rm~\ 
ivi_l/ivi \,±o\J cui 1 


R 


La Silla, Chile 


GROND (220 cm) 


g, r, i, z, J, H , K 


L ' Ampolla^ 




R 


Lowell (Perkins)^ 


Perkins 


R 


T 1 1 1 1 VI r 1 ] o 1 Tiro VI TP 
1 J 1 J, .1, 1 II, ± til Weill 


^T,T (df\ rm\ 

UDl 1 1U t,lll J 




Roque, Canary Islands* 


KVA (35 cm) 


R, polarization (no filter) 


Roque, Canary Islands* 


LT (200 cm) 


R, H 


San Pedro Martir* 


(84 cm) 


R 


St. Petersburg, Russia* 


(40 cm) 


R 


Tijarafe* 


(35 cm) 


R 




Radio 




CARMA, USA 


(array) 


92.5, 227.5 GHz 


Mauna Kea, USA* 


SMA(8 x 6 m) 


230, 345 GHz 


Medicina, Italy* 


(32 m) 


5, 8, 22 GHz 


Metsahovi, Finland* 


(14m) 


37 GHz 


Noto, Italy* 


(32 m) 


43 GHz 


Owens Valley, USA 


OVRO (40 m) 


15 GHz 


UMRAO, USA* 


(26 m) 


5, 8, 14.5 GHz 


* GASP-WEBT 






Calar Alto data 


was acquired as part 


of the MAPCAT 


project http: //www. iaa.cs/~iagudo/rcscarch/MAP CAT 





photometry with a comparison star taken in the same 
frame of 3C 279. Its position is R.A.=12h56ml6.90s, 
Dec=-05°50'43.0"(J2000) and its m agnitudes are V = 
13.660, J = 12.37 7 and Ks =11.974 (IRaiteri et al.lU99H: 
iCutri et al.ll2003l) . The photometric data have been cor- 
rected for the Galactic extinction with Ay = 0.093, 
Aj = 0.026, and A Ks = 0.010. The measured op- 
tical and near-IR magnitudes by Kanata in the V, J 
and Ks bands during the 2-year observations spanned 
m v = 15.54-17.27 (56 data points), mj = 13.00-14.58 
(37 data points) and m Ks = 11.21 - 11.47 (17 data 
points) . Those data points are also shown in Figure [7J 

We confirmed that the instrumental polarization was 
smaller than 0.1% in the Fband using the observations 
of unpolarized standard stars. Hence, we did not 
apply any corrections for it. The zero point of the 
polarization angle is corrected as standard system (mea- 
sured from north to east) by observing the po larized 
stars, HD 19820 and HD25443 (| Wolff et all Il99l . The 
polarization shows clear variability and the degree of 
polarization was measured in the range of 3 - 36 % dur- 
ing our 2-year observational campaign. As we reported 
in Paper I, we found a rotation of the polarization angle 



by 208° together with a sharp drop of the degree of 
polarization from ~ 30% down to a few %. The event 
was coincident with a 7-ray flare (Period D-E). In the 
second half of the 2-year observations, the source was 
generally in a quiet state in the optical band, and the 
degree of polarization was also relatively low. 



3.8. Optical and Near-Infrared: GROND 

The Gamma-R ay burst Optical/N ear-Infrared De- 
tector (GROND; iGreiner et all 12008ft mounted at the 
MPI/ESO 2.2 m telescope at LaSilla observatory in Chile 
observed the field of 3C 279 in two nights of 2008 July 
(2008 July 30 and 2008 July 31) and four nights in 2009 
January (2009 January 19 to 2009 January 22). In each 
observation, a total of 4 images in each g'r'i'z' filter with 
integrations times of 35 s and 24 images of 10 s exposure 
in each JHK S were obtained simultaneously. 

GROND optical and near-IR data were redu ced in 
standard manner using pvraf/IR AF (iTodvl 11993ft sim i- 
lar to the procedure outlined in iKriihler et al.l ( 2008). 
The stacked images of each observation were flux cal- 
ibrated against GROND observations of SDSS fields 



12 



Hayashida et al. 



(Abazaiia n et al.ll2009f ) taken immediately before or af- 
ter the field of 3C 279 for the optical g'r'i'z', and m ag- 
nitudes of 2MASS field stars (jSkrutskie et al.ll2006D for 
the JHK S filters. All data were corrected for the ex- 
pected Galactic foreground reddening of Er B _ V \ = 0.029 
according to iSchlegel et al.1 (|1998f ). Results of GROND 
observations are summarized in Table El 

3.9. Infrared: Spitzer 

We observed 3C 279 with Spitzer Infrared Spectro- 
graph (IRS), Multiband Imaging Photometer for Spitzer 
(MIPS) and Infrared Array Camera (IRAC) at several 
epochs in 2008 and 2009 under the Spitzer program 
PID50231 (PI A. Wehrle; see TableEJ). The observations 
were conducted once with each instrument in 2008 July 
and August and approximately daily during the instru- 
ment campaigns in the 2009 February-March visibility 
window. 

For the IRS observations, high-accuracy blue peakup 
observations on a nearby star were used to center the 
spectrograph slit on the target. 3C 279 was observed 
with the low-resolution SL2, SL1, LL2, and LL1 modules, 
for 3 cycles of 14 seconds at each of two nod positions. 
Spitzer-IRS data reductions began with S18.7 Spitzer Sci- 
ence Center pipeline-processed, background subtracted 
data. The background was removed by subtracting the 
alternate nod for each pointing. Additional processing 
steps were applied to clean bad data, remove fringes, and 
match and trim spectral orders. First, we cleaned bad, 
rogue pixels using the Spitzer Science Center procedure 
IRSCLEAN V2.0. One-dimensional spectra were then 
extracted using the standard point-source aperture and 
flux calibration in SPICE ver. 2.3. We used a custom 
spectral defringing tool to remove fringes introduced by 
the pointing-dependent instrumental flatfield. This tool 
uses a predetermined flatfield fringing correction func- 
tion, which is shifted to match and remove the observed 
fringes in the spectrum. Spectral orders were trimmed, 
and the SL2 and SL1 orders were scaled up by a factor of 
1.06 to empirically correct for pointing-dependent point- 
source slit losses. Finally, the nod-spectra were averaged 
and combined into a single spectrum covering 5.2-35 /im 
rest wavelength. Figure [3] shows the reduced IRS spectra 
in the vF v representation. 

We used the pipeline MIPS images (ver. 18) for aper- 
ture photometry using 13", 35" and 50" radius for 24, 70 
and 160 /im bands, respectively, with aperture correc- 
tions from Tables 3.13, 3.14 and 3.16 of the MIPS Data 
Handbook (ver. 3.2). No 160 /im data were obtained 
in August 2008 because the array was not cooled dur- 
ing that campaign. 3C 279 has very low ecliptic latitude 
(0.2°), hence, the observed transients can be attributed 
to passing asteroids which appeared in various MIPS im- 
ages. We used 6" radius apertures for IRAC photome- 
try on the pipeline data (ver. 18) with aperture correc- 
tions tabulated in Table 5.7 of the IRAC Data Handbook 
(ver. 3.0). 

The Spitzer-MlPS photometric repeatability and abso- 
lute calibration uncertainties at 24 /im are respectively, 
0.4% and 4%; a t 70 /im, 4.5 % and 5%; a nd at 160 um 
5 % a nd 12 % (lEngelbracht et al.l [20071: iGordon et all 
120071 iStansberrv et al.ll2007D . We therefore adopt over- 
all uncertainties of 10 %, 10 % and 20 % at 24, 70 and 
160 /im, respectively. No color correction has been ap- 
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Figure 9. Spectral energy distribution of 3C 279 in the infrared 
band measured by Spitzer-IRS. IRS spectra from highest to lowest 
are on 2009 March 3, 2009 March 4, 2009 March 5, 2009 March 6, 
2009 March 7, 2009 March 9, and 2008 August 16. Representative 
flux data points measured by Spitzer-MIPS (circles) and Spitzer- 
IRAC (triangles) are also included. 

plied to the data because the slope and smoothness of the 
spectrum over the bandpasses are not known. Figure [TU] 
describes the MIPS flux history during the 6 epochs from 
15 to 20 February 2009 together with i?-band flux for 
comparison. No significant flux variation is found in all 
MIPS bands during those epochs, which include period 
D. 

The Spitzer-IRAC calibration uncertainty is 3 % overall 
and h as photometric repeatability of 1.5% ()Reach et al.l 

2005) . We adopt the overall IRAC calibration uncer- 
tainty of 3%, but note the following characteristics of 
our images. In our IRAC frames, two standard compar- 
ison stars used in blazar monitoring were visible in the 
3.6 /im images (Star 1 and Star 2)0 One comparison 
star, Star 2, was visible in the 4.5, 5.8, and 8 /im images, 
located at the interstice of the chopping regions where 
the data are noisier than elsewhere. The spacecraft ori- 
entation, and hence the chopping orientation, was 180 
degrees different between 2008 July-August and 2009 
March. The standard deviations in comparison Star 2's 
measurements in 2009 March at 3.6, 4.5, 5.8 and 8 /im 
are 0.68 mJy, 0.08 mJy, 0.08 mJy and 0.09 mJy (5 %, 1 %, 
2 % and 3 %), respectively. The high 3.6 /im standard de- 
viation was affected by a single high value on 11 March 
2009, for which we found no obvious cause; excluding 
that value resulted in a standard deviation of 0.06 mJy 
(0.5%). In contrast, the flux of 3C 279 shows a steady 
decrease of 10%, 12%, 14% and 13% at 3.6, 4.5, 5.8 
and 8.0 /im, respectively, during the 6 epochs from 2009 
March 10-16 as shown in Figure [TT1 

3.10. Radio: CARMA 

Observations were obtained at mean frequencies of 92.5 
and 227.5 GHz using the Combined Array f or Research 
in Millimeter-wave Astronomy (CARMA; Boc k et all 

2006) . In all cases the nominal signal to noise ratio 
exceeded 400 and calibration uncertainties dominated 

65 sec http : //www . lsw .uni-heidelberg . de/projects/extragalactic/charts, 
IRaiteri et af.l UWb]). IVillata et al.1 dlWTI) and 

http : //quasar . Colgate . edu/~ tbalonek/optical/3C279compstars . gif 
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Table 3 

Results of GROND observations. 





2008 July 31 


2009 January 19-22 




filter 


AB magnitude 


flux [mjy] 


AB magnitude 


flux [mjyj 


A/A{V) a 


g' 


17.62 ± 0.05 


0.324 ±0.015 


16.06 ± 0.05 


1.37 ±0.06 


1.23 


r> 


17.24 ±0.05 


0.462 ± 0.021 


15.49 ±0.05 


2.31 ±0.11 


0.80 


i> 


16.83 ± 0.05 


0.671 ± 0.031 


15.05 ± 0.05 


3.47 ±0.16 


0.62 


z' 


16.67 ±0.05 


0.776 ± 0.036 


14.80 ±0.05 


4.37 ±0.20 


0.45 


J 


16.05 ± 0.06 


1.387 ±0.077 


14.16 ± 0.06 


7.86 ± 0.44 


0.29 


H 


15.60 ±0.07 


2.098 ±0.133 


13.62 ± 0.07 


12.90 ±0.81 


0.18 


K 


15.18 ± 0.09 


3.062 ± 0.249 


13.19 ± 0.09 


19.18 ± 1.56 


0.14 



Note. — Results of both AB magnitude and flux are corrected for Galactic extinction. 
No significant daily variability was observed during the observations in 2009 January. 
a dereddening factors for correction of Galactic extinction. 



Table 4 

Spitzer observation 



Instrument 


Start time (UTC) 


(MJD) 


Duration (min) 


ObsID 


MIPS 


2008.7.31 11:09:06.3 


54678.4647 


9.19 


27434240 




2009.2.15 07:08:08.6 


54877.2973 


14.61 


27438592 




2009.2.16 23:12:32.0 


54878.9670 


14.61 


27438080 




2009.2.17 20:42:36.5 


54879.8629 


14.61 


27438848 




2009.2.18 19:39:30.8 


54880.8191 


14.60 


27438336 




2009.2.19 14:19:04.4 


54881.5966 


14.60 


27439360 




2009.2.20 09:23:37.4 


54882.3914 


14.60 


27439104 


IRS 


2008.8.16 14:32:42.6 


54694.6060 


17.92 


27425024 




2009.3.3 12:22:35.7 


54893.5157 


18.15 


27435776 




2009.3.4 19:53:03.5 


54894.8285 


18.15 


27437312 




2009.3.6 10:50:02.7 


54896.4514 


18.16 


27436544 




2009.3.7 00:46:03.4 


54897.0320 


18.16 


27435520 




2009.3.8 12:08:53.5 


54898.5062 


18.18 


27437056 




2009.3.9 18:39:51.2 


54899.7777 


18.18 


27436288 


IRAC 


2008.8.17 06:57:26.3 


54695.2899 


10.64 


27429632 




2009.3.10 19:26:05.5 


54900.8097 


10.69 


27433216 




2009.3.11 22:05:54.7 


54901.9208 


10.70 


27432448 




2009.3.13 05:13:56.9 


54903.2180 


10.71 


27433728 




2009.3.14 04:12:41.8 


54904.1755 


10.72 


27432960 




2009.3.15 03:04:16.9 


54905.1280 


10.73 


27433984 




2009.3.16 14:58:46.6 


54906.6241 


10.75 


27433472 
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Figure 10. Light curve of 3C 279 at 24 /an, 70 fim and 160 fim 
measured by Spitzer-MIPS. The error bars correspond to 10%, 10% 
and 20 %, re spectively, for each band as mentioned in the text in 
Section 13.91 Optical ij-band data taken by the ground- based tele- 
scopes are also plotted in red color for comparison. 



the errors. The source is bright enough to permit self- 
calibration on timescales of less than a minute and so 
atmospheric decorrelation was not expected to affect our 
results significantly even at the long baselines. However, 



Figure 11. Light curve of 3C 279 at 3.6 fim, 4.5 (im, 5.8 fim and 
8.0 fim measured by Spitzer-IRAC. The error bars correspon d to 
3% for all IRAC bands as mentioned in the text in Section 13.91 
Optical i?-band data taken by the ground-based telescopes are also 
plotted in red color for comparison. 

observations in poor weather were not used due to the 
difficulty of reliably measuring pointing offsets in these 
conditions. 

Data calibration and a nalysis was don e with the 
MIRIAD software package (Sau lt et allll995f ). Flux den- 
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sities were determined by first using phase self calibra- 
tion with a short enough averaging interval to avoid any 
atmospheric phase de-correlation, then the flux density 
was determined from the vector average fringe ampli- 
tude at the position of 3C 279 over all baselines. For 
a strong point source such as 3C 279 this provides very 
robust and unbiased amplitude estimate independent of 
the weather or the interferometer baselines. We rely 
on regular system temperature measurements to provide 
flux calibration relative to the the fixed system sensi- 
tivity. The absolute flux calibration of CARMA obser- 
vations is usually quoted as 10-15%. However, based 
on measurements made on the blazar 3C 454.3 we esti- 
mated the relative flux calibration at each frequency to 
be within 5 % at 3 mm and 10 % at 1 mm. The radio 
fluxes at 92.5 and 227.5 GHz measured by CARMA cor- 
respond to i*92.5 GHz = 11-7 — 19.7 Jy (14 data points) 
and -F227.5GHZ = 6.3 — 9.2 Jy (14 data points). Figure |H] 
includes the flux history of those radio data. 

3.11. Radio: OVRO 40m 

The Owens Valley Radio Observatory (OVRO) 40 m 
radio data were collected as part of an ongoing long- 
term, fast-cadence 7- ray blazar monitorin g campaign, 
described in detail in iRichards et al.l (|2011f) . Flux den- 
sities were measured in a 3 GHz bandwidth centered on 
15.0 GHz using dual, off-axis 2'.5 FWHM beams with 
12f95 separation. Dicke switching against a blank sky 
reference field to remove gain fluctuations and atmo- 
spheric and ground contamination were used. Flux den- 
sities from this program are found to have a minimum 
uncertainty of 4 mJy (mostly thermal) and a typical un- 
certainty of 3 % for brighter sources. During the period 
included here, 3C 279 was observed as a pointing cali- 
brator. The flux density scale was referred to th e value 
for 3C 286 (3.44 Jy at 15 GHz: lBaars eTa^[l977t) with a 
scale uncertainty of about 5%. The radio flux at 15 GHz 
measured by OVRO was ranging from 11.1 to 18.0 Jy 
among 124 data points during the 2-year observations. 
The light curve of the OVRO radio data is also plotted 
in Figure HI 

4. RESULTS OF THE MULTI- WAVELENGTH 
OBSERVATIONS 

4.1. Correlations of light curves in various bands 

The multi-band light curves of 3C 279 are presented in 
Figure [T3] They include (a) 7-ray flux above 200 MeV 
(Fermi-L AT), (b) 2-10 keV X-ray flux measured by Swift- 
XRT and RXTE-PCA, (c) optical-UV fluxes in i?-band 
(GASP), F-band (Si/ra/t-UVOT and Kanata) and W2- 
band (Swift-WOT), (d,e) degree and angle of optical 
polarization (Kanata and KVA) and (f) radio fluxes in 
the 230, 37, 15 and 5 GHz bands (GASP, CARMA, and 
OVRO). We note that the X-ray fluxes determined by 
Suzaku and XNINI-Newton are entirely consistent with 
those plotted in Figure []~~] The extensive data set ob- 
tained in many bands for 3C 279 allows us to make gen- 
eral statements regarding the relative flux variability in 
various spectral bands, and the relationship of the time 
series to each other. The first such feature of the multi- 
band light curves is a general - although not exact - trend 
where the IR through optical emission seems to be cor- 
related with the 7-ray flux. We calculated the Discr ete 
Correlation Function (DCF; lEdelson fc Krolikl [l988h to 



quantify the correlation of the flux variations between 
the 7-ray and other bands, and to determine whether we 
can measure any clear lag between the bands. 

For the DCF calculations, we use the 7-ray fluxes aver- 
aged over an interval of 1 day as shown in the top panel of 
Figure [T] The resulting DCF between 7-ray and optical 
i?-band fluxes is shown in Figure [T3] Positive values of 
V correspond to flux variations in the 7-ray band lagging 
flux variations in the other bands. In the DCF between 
7-ray and optical i?-band fluxes, a peak can be seen close 
to zero lag. We fit the DCF data points in the range be- 
tween —30 and 5 days using a Gaussian function of the 
form DCF(t) = C max x exp[(r - r ) 2 /cr 2 ], where C max is 
the peak value of the DCF, tq is the time at which the 
DCF peaks, and a is the Gaussian width of the DCF. 
The fit yields a position of the peak at tq = —10.7 ± 0.7 
days, corresponding to a value of C max = 1.07±0.03 with 
a dispersion of a = 19.4 ±1.4 days. The result implies 
that the optical emission is possibly delayed with respect 
to the 7-ray emission by about 10 days. 

In the framework of the one-zone synchrotron + 
external-radiation Compton (ERC) models, the same 
electron population, of roughly the same energies, is re- 
sponsible for the radiation in both the optical and 7-ray 
bands. There, the observed lag can result from different 
profiles of the decreasing magnetic and radiation energy 
densities along the jet: we show that idea quantitatively 
in Appendix A. As is shown there, a very steep drop of 
the external radiation energy density is required to ex- 
plain the lag in a conical jet with magnetic field B' cx 1/r 
where r is the distance along the jet. This condition can 
be relaxed i n the scenario involving the re- confinement 
of a j e t fe.g..lDalv fc Marscherlll9 88; Ko missarov fc Fallel 
119971: iNalewajko fc Sikoral I2009Q . In such a case, the 
magnetic field intensity can drop more slowly than 1/r. 
If the lag of the optical emission is confirmed, the ap- 
plication of the results in Appendix A to the ~ 10-day 
lag may imply the location of the active "blazar zone" 
at distances of a few pc in agreement with those pos- 
tulated to explain the optical polarization swing (Pe- 
riod D-E) in terms of a region containing an enhanced 
density of ultra-relativistic electrons propagating along a 
curved trajectory (Paper I). It is worth noting that simi- 
lar 7-ray/optical lags have be en reported during the o ut- 
bursts of 3C 279 in earl y 1999lHartman et all (|2001bl) . of 
PKS 1502+106 in 2008 (lAbdo et al.l2010ch.~of PKS 1510- 
089 i n early 2009 (|Abdo et al.ll2010etlD'Ammando et al. 
20T]ll and of AO 0235+164 i n late 2008 (|Agudo et al. 
2011al lAckermann et aLll2012f) . On the other hand, no 
significant lags between 7-ray and optic al signals have 
been detected in 3C 454. 3 in late 2008 (Bonnin g et al.l 
I2009t Uorstad et all 12010ft . in 3C 66A in 2008 Octo- 
ber (lAbdo et all 120111) a nd in OJ 287 in 2009 Octo- 
ber (jAgudo et al.ll2011b[ ). Based on investigations of 
long-term light cur ves of 3C 454.3 during 2008-2010, 
iRaiteri et al.l ([201 If ) have shown that the optical and 7- 
ray flux variations are not always simultaneous and have 
proposed a geometrical scenario to explain the change 
in the 7/optical flux ratio during the outburst peaks in 
3C 454.3. It is expected that the on-going multi-band 
monitoring of blazars will enable us to quantify such lags 
and find out how common they are. 

Different behavior is apparent in the radio flux, where 
the energies of radio-emitting electrons are very different 
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Figure 12. Multi-band light curves of 3C 279 for 2 years from 2008 August to 2010 August, (a): Gamma-ray flux above 200 MeV 
averaged over 3 days, (b): X-ray flux between 2 and 10 keV measured by Suii/i-XRT and RXTE-PCA. (c): UV-optical fluxes in i?-band 
(red), V-band (blue) and 1^2-band (magenta), (d): Polarization degree in the optical band, (e): Polarization angle in the optical band. 
The horizontal dashed lines refer to the angle of 50° and —130°. (e): Radio fluxes in 230 GHz band (magenta), 37 GHz band (orange), 15 
GHz band (blue) and 5 GHz band (green). All X-ray, UV and optical data are corrected for the Galactic absorption. 

from the energies of the electrons involved in produc- 
ing the observed optical and 7-ray emission. Variabil- 
ity appears to be much less rapid, and the excess vari- 
ance (-Fvai-; see definition in Eq. [1} in the radio regime 
is quite modest: for instance, 0.145 ± 0.004 at 37 GHz, 
0.165 ± 0.001 at 15 GHz and 0.104 ± 0.001 at 5 GHz. 
Those values are significantly less than ones in the 7- 
ray or optical bands. This suggests that the synchrotron 
emission from the 7-ray emitting region is self-absorbed 
at these wavelengths. The observed radiation is produced 
at much larger distances, where the light-travel effects 
smear out the sharp, rapid variability patterns observed 
in the optical and 7-ray bands. 

Perhaps the most surprising behavior - and difficult to 
explain in the context of simple, one-component, single- 
zone models - is the relationship of the X-ray light curve 
to those in the IR-optical or 7-ray bands. In Paper I, we 
reported that the X-ray time series exhibits a relatively 
rapid, symmetrical flare at ~ MJD 54950 (Period F) with 
a duration of ~ 20 days, which is not accompanied with 
any prominent IR/optical or 7-ray flares. As we argued 
in Paper I, the hard (rising in vF v representation) X-ray 
spectrum is unlikely to be the "tail" of the synchrotron 
emission, but instead, it is more likely to be produced by 
the low-energy end of the electron distribution radiating 
via inverse Compton process. 

The continuing monitoring of the object in the X-ray 
band revealed another X-ray flare at ~MJD 55040 (Pe- 
riod G), ~ 90 days after the first X-ray flare. The sep- 
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Figure 13. Discrete Correlation Function (DCF) derived for the 
7-ray and optical _R bands. Positive values of V correspond to flux 
variations in the 7-ray band lagging flux variation in the optical 
band. The red curve represents a Gaussian fit to the data between 
—30 and 5 days. See the text for the fitting results. 
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Figure 14. DCF between derived for the 7-ray and X-ray bands. 
Positive values of V correspond to flux variations in the 7-ray band 
lagging flux variation in the X-ray band. 
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aration of the two X-ray flares is remarkably close to 
the temporal separation of the two 7-ray flares, with the 
two pairs delayed with respect to each other by ~ 155 
days. Figure [14] presents the calculated DCF between 
7-ray and X-ray fluxes, which shows a modest peak at 
~ —155 day with a correlation coefficient of 0.6-0.7 and 
indicates no correlation between the the 7-ray and the 
X-ray bands with zero lag. While confirming the phys- 
ical connection of the two pairs would be very impor- 
tant, we cannot currently envision any situation where 
the two would be causally connected: the 155-day lag 
would imply the distance of the X-ray flare production 
~ 155r? light days ~ 50(rj/20) 2 pc and at such a dis- 
tance should be accompanied by radio flares, which are 
not seen in our data. In such a scenario, the X-ray flares 
should be significantly broadened compared to the 7-ray 
flares, however, we observe a similar temporal structure 
in both bands. Furthermore, we note that there are some 
optical and 7-ray peaks that might well be associated 
with the second X-ray flare. Hence, it is possible that 
the two prominent 7-ray/optical flares (Period B and D), 
together with the subsequent two X-ray flares (Period 
F and G), form a sequence of 4 events separated by a 
similar time intervals. Those intervals, in turn, can be 
possibly determined by instabilities in the jet launching 
region. Here, the different broad band spectra during 
these events may result from small changes of param- 
eters, such as the jet direction, Lorentz factor, and/or 
location and geometry of the dissipation event. 

A weak (and sporadically almost absent) correlation 
between X-rays and other spectral bands can also result 
from such processes which preferably contribute to ra- 
diation in the X-ray band. They can be related to the 
following three mechanisms/scenarios: 

1. Bulk-Compton process. This involves Compton- 
scattering of ambient optical/UV light by the cold 
(non-relativistic) electrons in the jet. This mecha- 
nism is most efficient close to the accreting black 
hole where the processes responsible for the vari- 
ability of X-rays may operate independently of 
those at larger distances and p roducing there vari- 
able non-thermal radiation (Bcgclman fc Sikoral 
1987). A drawback of this scenario can be that 
the bulk-Compton spectrum is predicted to have a 
similar shap e as the spectrum of th e external ra- 
diation field (jAckermann et aLll2012f) , which signif- 
icantly differs from what we observe in the X-ray 
band. 

2. Inefficient electron acceleration. Acceleration of 
the relativistic electrons at proton-mediated shocks 
is likely to proceed in two steps: in the first one low- 
energy electrons may be pre-accelerated via, for ex- 
ample, some collective processes involving protons; 
in the second step, they may participate in the first- 
order Fermi acceleration process. If under some 
conditions the electron-proton coupling is ineffi- 
cient, the fraction of electrons reaching the Fermi 
phase of acceleration will be small. In this case the 
X-rays, originating from lower energy electrons, are 
produced efficiently, while the 7-rays and optical 
radiation which involve more relativistic electrons 
- are not. 



3. The X-rays can be also contributed by hadronic 
processes, specifically by the pair cascades powered 
by proto ns losing their energy in the pho to-mesonic 
process ([Mannheim fe; Biermannl Il992f ). For this 
proce ss to be efficient, it requires e xtreme condi- 
tions ([Sikora e t al. 2009; Si koral2011[ ). however, op- 
erating in the very compact central region, at dis- 
tances less than few hundred gravitational radii, it 
may occasionally dominate in the X-ray band. 

4.2. Broad-band Spectral Energy Distribution 

Figures p] and M show broad-band SEDs of 3C 279 
in all periods as defined in Table [1] In addition, wc 
also extracted a SED using data taken on 2008 July 
31 (MJD 54678), which has a good energy coverage of 
the synchrotron emission component including Spitzer 
and GROND data, although the 7-ray data by Fermi- 
LAT are not available at that time because this was be- 
fore the beginning of normal, all-sky science observations 
with Fermi-L AT. Both SEDs for Period D (2009 Febru- 
ary: corresponding to the brightest 7-ray flare coinci- 
dent with the optical polarization swing) and Period F 
(2009 April: corresponding to the first isolated X-ray 
flare) have already been partially reported in Figure 2 
of Paper I. New Spitzer-MIPS data points are included 
in the SED for Period D in this paper. In Period C, 
there are observations by MAGIC, which provide upper 
limits above 100 GeV ([Aleksic et al.ll2011[ ). For compar- 
ison, we also include very-high-energy 7-ray fluxes de- 
tected with the MAGIC telesc ope in 2006 February as 
gray points (|Albert et al.ll2009[ ) in Figure [15l 

This is the richest set of time-resolved spectra ever 
collected for this source. The spectral coverage of the 
synchrotron bump is unprecedented, allowing us not only 
to constrain the parameters of the emission models, but 
also to study their time evolution. As we discussed in 
Sect ion |2~41 the shape of the 7-ray spectrum deviate from 
a simple power law, in similarity to other FSRQ blazars. 
Strong variability, over one order of magnitude, is evident 
in nearlR/optical/UV and 7-ray bands. This contrasts 
with the moderate variability in the radio/mm and X-ray 
bands. 

Particularly interesting is the behavior of this source 
in the mid-IR band, around ~ 10 13 Hz, where signif- 
icant spectral variability is observed. In the low state 
in Period A the mid-IR spectrum is relatively soft and 
can be extended with a power law shape to the opti- 
cal/UV band. In this case, the synchrotron component 
peaks in the mm/sub-mm band (~ 10 11 — 10 12 Hz). How- 
ever, in the high state in Periods D and E, the mid-IR 
spectrum is much harder and shows a significant cur- 
vature. In Period D, there is a clear spectral break at 
~3x 10 12 Hz (~ 100 jitm). The spectral index between 
the 70 fim (~ 4.3 x 10 12 Hz) and 160 fim (~ 1.9 X 10 12 Hz) 
points is «7o-i6o = 0-35 ± 0.23, taking into account sys- 
tematic errors described in Section [3.9l The synchrotron 
peak is located in the mid-IR band, at a frequency one or- 
der of magnitude higher than in the low state. This indi- 
cates that there are two independent synchrotron bumps, 
possibly produced at different locations. The mid-IR- 
peaking component, seen only in the IR/optical/UV flar- 
ing state, is characterized by a strong and rapid variabil- 
ity. The mm/sub-mm peaking component is more persis- 
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Figure 15. Time-resolved broad-band spectral energy distributions of 3C 279 measured in Period A-H (as denned in Table[JJ and on 2008 
July 31 (MJD 54678), covered by our observational campaigns in 2008-2010. X-ray, UV-optical-nearlR data are corrected for the Galactic 
absorption. 5-digit numbers in the panel indicate MJD of the periods. For comparison, th e gray open circles in the very-high-energy 7-ray 
band represent measured spectral points by MAGIC in 2006 February (Albert ct al. 2003). 



tent and dominates when the source is in the low state. 
The complex shape of the SED between the mm band 
and the 70 /im point in Period D requires a coexistence 
of these two components. A similar scenario of multiple 
synchrotron component s was investigated in the case of 
3C 454.3 bv lQgfe et all (poll . 

In the X-ray band, despite the smaller variability am- 
plitude, we observe some spectral changes. In particular, 
in Periods F and G, which represent the two isolated X- 
ray flares, the spectrum is very similar and harder than 
on average. Figure [15] shows that these flares are not en- 
ergetically important. If we extrapolate the X-ray spec- 
tra with power laws to the 7-ray band, we under-predict 
the observed 7-ray flux in Periods B, C, D and H. Periods 
B, C, D coincide with the high-activity 7-ray state. This 
indicates that the X-ray flux cannot originate from the 
same emission component as the 7-ray flux, at least in 
the flaring state. Because the 7 rays are correlated with 
the optical flux but not with the X-ray flux, the 7 rays 
can be related to the mid-IR-peaking synchrotron bump 
while the X-rays may correspond to the mm/sub- mm 
peaking synchrotron bump. We explore this possibility 
when modeling the SEDs at Periods A and D in Section 

5. MODELING THE BROAD-BAND EMISSION 

We have fitted sele cted SEDs with one-zo ne leptonic 
models described in iModerski et al.l (|2003l ). including 
synchrotron emission and self-absorption, Comptoniza- 
tion of the local synchrotron radiation [SSC component] 
and external photons [ERC component] , but also includ- 
ing the opacity due to internal pair-production. The ex- 



ternal radiation includes broad emission lines (BEL) and 
infrared dust emission (IR). Their energy densities in the 
jet co-moving frame as functions of the distance r along 
the jet are approximated by the formulae: 



*BEL 



(r) 



,BEL 



3t"b EL c[1 + {r/rBEh) PBEL } 



Z*r 2 m c[l + (r/nR)fl«] 



(2) 



(3) 



where £bel = 0.1 and £ir = 0.1 are the fractions of the 
disk luminosity i D ~ 2 x 10 45 ergs -1 reprocessed into 
emission lines and into hot dust radiation, respectively, 
t-bel = 0.1(Lzx46) 1/2 pc and r m = 2.5(L DA6 ) 1 / 2 pc 
[£d,46 = ^d/10 46 ] are the characteristic distances where 
such reprocessing takes place, and Tj is the jet Lorentz 
factor. The external radiation fields are approximated in 
the jet co-moving frame by Maxwellian spectra peaked 
at photon energies -E BEL ~ 10 eV x Tj and ~ 
0.3 eV x Ij. While the radiation density profile in the 
frame external to the jet should satisfy /3bel(ir) < 2, it 
is not applied to the profile in the jet co-moving frame. 
This is because the transformation of radiation density 
depends on the angular distribution of external radiation, 
with radiation arriving at small incidence angles to the 
jet velocity vector being actually deboosted. This can re- 
sult in a steeper profile of the radiation densit y in the jet 
co-m oving frame. Here, we ado pt /?b el = 3 (jSikora et al.l 
2009) and /3tr = 4 (see Section [5TTj) . We assume a coni- 
cal jet geometry; the magnetic field, assumed to be domi- 
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Figure 16. Time-resolved broad- band spectral energy distributions of 3C 279 in each period (A-H), covered by the campaigns. The data 
points are the same as ones in Fig. 1151 but are plotted in a separate panel for each period. 5-digit numbers in each panel indicate M JD of 
the observing period of each broad-band spectrum. 
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nated by the toroidal component, is taken to decline with 
distance r as B' oc 1/r. Electrons are injected with a 
double-broken power law distribution Q(j) oc 7~ p with 
p = pi for 7 < 7b r i, p = p 2 for 7 br i < 7 < 7 br 2 and 
p = P3 for 7 > 7br2- Their evolution, including injec- 
tion at a constant rate as well as radiative and adiabatic 
cooling, is followed over a distance Ar = r/2, where r is 
the position at which the injection ends. The emission is 
integrated over spherical thin shells within a conical re- 
gion of opening angle 0j = 1 /Tj . The observer is located 
within the jet opening cone, i.e. (9 b s < 0j- 

We begin by modeling the SED in Period D, which is 
the highest 7-ray state reached by the source during our 
observational campaigns. In Paper I, we showed that 
the flare event was accompanied by an optical polariza- 
tion swing and proposed two interpretations of this event. 
The first one involved a cloud containing ultra-relativistic 
particles propagating along a curved trajectory. The du- 
ration of the polarization swing constrains the location 
of the cloud to be at a few parsecs from the central su- 
permassive black hole, in the region where external ra- 
diation is dominated by the infrared dust emission. In 
Section I5.1[ we present an ERC-IR model describing the 
SED in Period D and a physically related model of the 
SED in Period E. The second interpretation of the polar- 
ization swing involved the jet precession, which allowed 
arbitrary location of the emitting region, including the 
broad-line region. In Section 15.21 we present an ERC- 
BEL model of the SED in Period D. We show that in 
this scenario the far-IR break arises due to synchrotron 
self-absorption. We also show an ERC-IR model of the 
SED in Period A, which can explain the mm/far-IR and 
X-ray emission, as well as the low-state optical and 7-ray 
flux levels. 

We assume the scenario where the X-ray emission is 
unrelated to the flaring component, since it showed lit- 
tle variability during the correlated 7-ray/optical flares. 
Our one-zone models of the flaring states are fitted only 
to the IR/optical/UV and 7-ray data, treating the simul- 
taneous X-ray spectrum as only an upper limit to the 
SSC component and the ERC component from the low- 
energy electrons. The large 7-ray/X-ray luminosity ratio 
forces us to adopt a very hard electron energy distribu- 
tion at low energies (p — 1), which can be alternatively 
obtained by imposing a minimum electron Lorentz factor 
7min > 1- 

5.1. Propagation scenario for the emitting region 

An intrinsically spherically symmetric emitting region 
is expected to produce the observed electric polarization 
vector aligned with the p rojected velocity of the emitting 
region. Nalcwaiko (2010) presented a simple model of its 
trajectory to explain the event of simultaneous smooth 
variations of the polarization degree and angle during the 
polarization swing which has been reported in Paper I. 
This model adopts a constant jet Lorentz factor Tj = 15 
and can be used to predict the viewing angles for a given 
observation time. For Period D we estimate 6 bs.B ~ 
1.5°, while for Period E: 6> bs,E ~ 2.4°. Between Periods 
D and E (At ~ 15 days) the emitting region propagates 
over a distance Ar ~ r?c(Ai) ~ 2.8 pc. 

In Figure [T7] we show Model Dl fitted to the SED in 
Period D at r = r m and Model El fitted to the SED in 



Period E at r = riR+Ar. Model parameters are listed in 
Table [5l Both models use the magnetic field scaled to the 
same value at the distance of 1 pc. In order to explain the 
difference in the luminosity ratio of the ERC component 
and the synchrotron component, which decreased by fac- 
tor ~ 4 between Model Dl and Model El, we assume a 
distribution of the co-moving IR radiation energy den- 
sity dropping steeply with distance, adopting flm ~ 4. 
This corresponds to a strongly stratified torus structure, 
with a significant concentration of hot dust very close t o 
the sublimation radius (see, e.g.. IMor fc Netzerl 12012ft . 
We should note that, although the relatively soft 7-ray 
spectrum was observed at Period E (r = 2.64 ± 0.32, 
see Table [lj , the peak of the ERC-IR component in the 
Model El falls at - 800 MeV, in the Fermi-LAT band. 

The far-IR spectral break in Period D requires a 
sharp break in the electron distribution function at 
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Figure 17. Emission models Dl (red line) and El (blue line) fitted 
to the spectral states at Periods D and E, respectively. Periods D 
and E correspond to the first 5 days and the last 3 days of the 
7-ray flaring event accompanied by an optical polarization change, 
respectively. Those models adopt our 'propagation scenario', where 
external radiation is dominated by infrared dust emission. We 
assume the X-ray emission is not related to the flaring component, 
and consider the X-ray fluxes as only upper limits to the SSC and 
the ERC co mpon ents during the flaring event. See the text in 
Section l5l and 15711 for details of the models and Table [5] for model 
parameters. 
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Figure 18. Emission Models A2 (green line) and D2 (red line) fit- 
ted to the spectral states at Periods A and D, respectively. Period 
A represents a quiescent state, and Period D corresponds to the 
7-ray flaring event accompanied by an optical polarization change. 
Those models adopt our 'jet precession scenario', which assumes 
the 7-ray flaring event (Period D) occurs within the broad-line 
region while the low-steady emission component (Period A) is gen- 
erated outside the broad-line region. See the text in Section ISI and 
I5.2l for details of the models and Table [5] for model parameters. 
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Table 5 

Parameters of emission models. 



Model 


Dl 


El 


D2 


A2 


ext. rad." 


IR 


IR 


BEL 


IR 


r [pc] 


1.1 


3.9 


0.045 


3.9 


R [pc] b 


0.07 


0.26 


0.0023 


0.19 




15 


15 


20 


20 


flj h 


3.8 


3.8 


2.9 


2.9 


8 obs [°] 


1.5 


2.4 


1.7 


1.7 


B' lpc [G] c 


0.14 


0.14 


0.15 


0.15 


M cxt I 10 " 4 er § cm" 3 ] d 


78 


0.97 


8 x 10 4 


1.8 


7brl 


800 


800 


170 


440 


7br2 


6500 


5000 


1000 




Pi 


1 


1 


1 


2.2 


P2 


2.5 


2.6 


2.4 


3.4 


P3 


6 


4.2 


7 




a A dominant component 


of the 


external radiation. 


IR: 



infrared dust emission, BEL: broad emission lines. 

Radius of source emission region. 
c Magnetic field intensity at the distance of 1 pc. 

Energy density of the dominant component of the external 
radiation in the jet co-moving frame at the given distance 
r. 

Tbri = 800. As the cooling break is expected at j c ~ 
3m e c 2 / (2<jTRu' ext ) ~ 660 (where R is a radius of source 
emission region and u' ext is the energy density of the ex- 
ternal radiation in the jet co-moving frame), 7b r i is lo- 
cated just within the fast-cooling regime. The electron 
distribution in the fast-cooling regime cannot be harder 
than p = 2, hence the resulting synchrotron spectral in- 
dex a = (p — l)/2 should be larger than 0.5. In fact, the 
mean value of the observed spectral index between 70 /im 
and 160 /im is smaller than 0.5 (070-160 = 0.35 ± 0.23), 
which cannot be explained if the electron cooling is effi- 
cient. However, the uncertainty of the measurement does 
not allow us to reject this scenario. 

5.2. The jet precession scenario: two synchrotron 
components 

Alternatively, if the jet precession can cause the ob- 
served 7-ray flare event with the polarization swing, the 
7-ray/optical emission can be generated much closer to 
the central black hole, even within the broad line region 
(see also in Paper I). Therefore, we also attempted to 
model Period D placing the emitting region at Tbel ■ For 
Tj = 15, with model parameters fitted using the syn- 
chrotron and ERC components, the X-ray flux is over- 
produced by the SSC process. To alleviate this problem, 
we increased the jet Lorentz factor to Tj =20. In Figure 
[T8l we show Model D2 with parameters listed in Table El 
The magnetic field strength scaled to the distance of 1 pc 
is almost the same as the value in Model Dl. Because of 
a smaller size of emission region and higher energy den- 
sity of the locally produced synchrotron radiation, the 
synchrotron self-absorption is able to produce a spectral 
cut-off at a higher frequency of ~ 3 x 10 12 Hz (~ 100 ^m) , 
consistent with the far-IR break. This interpretation has 
an advantage that it also could explain the observed hard 
spectral index between 70 /mi and 160 /im, even smaller 
than 0.5, independently of details of the electron energy 
distribution. 

The low-energy synchrotron component, dominating 
the mm/sub-mm band, must be produced in a much 
larger region, placing it far outside the broad-line re- 
gion. In Figure [TH1 we present Model A2, fitted to the 
SED at Period A. We kept the Lorentz factor and the 



magnetic field consistent with Model D2, but we set the 
source at the distance ~ 4 pc, the same as in Model El. 
This low state model of Period A can reproduce both ob- 
served X-ray and 7-ray spectra by a single broken power 
law electron distribution. The 7-ray spectral index is 
consistent with the IR/optical/UV spectral index. The 
synchrotron self-absorption is effective at ~ 10 11 Hz and 
the spectral peak is located in the mm/sub-mm band. 

Those results suggest the existence of two synchrotron 
components: one peaking in the mm/sub-mm band and 
the other peaking in the mid-IR band. The component 
with the peak in the mid-IR band is more variable, and 
can be produced at shorter distances, within the broad- 
line region, where the far-IR break can be explained by 
the synchrotron self-absorption. 

6. CONCLUSIONS 

This paper reports details of the multi-band campaigns 
on the well-known blazar 3C 279 during the first two 
years of the Fermi mission between 2 008 and 2010. Some 
key results were already presented in lAbdo et al.1 ([2010a, 
Paper I). Most important of them was the coincidence 
of a dramatic 7-ray/optical flare with a change in the 
optical polarization, which we interpreted as the result 
of a compact emitting region: either propagating along 
a curved relativistic jet or located at a constant distance 
in a precessing jet. In addition, we reported on an "iso- 
lated" X-ray flare, an event without a clear counterpart 
in other bands, and taking place a few months after the 
7-ray/optical flare. Here, we extended the observation 
epoch until 2010 August yielding the best coverage of 
time-resolved SEDs ever collected for 3C 279 from radio 
through high-energy 7-ray bands. Based on those data 
we arrived at several new conclusions about the structure 
and emission models of the relativistic jet in the quasar: 

• In the high-energy 7-ray band measured by Fermi- 
LAT, the source exhibited two prominent flares 
reaching as high as ~ 3 x 10 -6 photons cm~ 2 s _1 
above 100 MeV in the first year while it was in a 
relatively quiescent state in the second year. No 
significant correlation between flux and photon in- 
dex has been measured in similarity to other LAT 
blazars. The 2-year averaged 7-ray spectral shape 
above 200 MeV clearly deviates from a single power 
law. The broken-power law model returns a break 
energy within a few GeV range, which does not ap- 
pear to vary with the source flux. Such behavior is 
similar to that observed in other bright FSRQs. 

• The superb temporal coverage allowed us to mea- 
sure in detail the cross correlation of the 7-ray 
and optical fluxes. The optical signals appear to 
be delayed with respect to the 7-ray signals by 
~ 10 days. Such a lag can be explained in terms 
of the simple synchrotron and inverse-Compton 
model, in the scenario where a cloud containing 
ultra-relativistic electrons propagates down the jet 
through the regions where the ratio of the external 
radiation energy to the magnetic energy densities 
decrease with distance. We have verified this idea 
qualitatively (see Appendix A), but it still needs 
specific numerical modeling to be confirmed quan- 
titatively. 



Radio to high-energy 7-ray observations of 3C 279 in 2008-2010 



21 



X-ray observations reveal a pair of pronounced 
flares separated by ~ 90 days. Those are not con- 
temporaneous with a pair of bright 7-ray /optical 
flares - also separated by ~ 90 days - but instead, 
are delayed with respect to the 7-ray/optical flares 
by about 155-days. Because of such a long delay, it 
seems implausible that these events are causally re- 
lated. Instead, the possible scenarios of the X-ray 
flares may involve changes of the source parameters 
such as the jet direction, Lorentz factor, and/or lo- 
cation of the dissipation event, or may require more 
'exotic' solutions, for instance: bulk-Compton pro- 
cess; inefficient electron acceleration above a given 
energy; hadronic processes. At this stage we can- 
not discriminate among any of those scenarios. 



The spectral coverage of the infrared band with 
Spitzer enabled us to probe the detailed structure 
of the low-energy spectral bump, attributed to the 
synchrotron radiation. Significant spectral vari- 
ability, with soft /power-law spectra in the low state 
and hard/curved spectra in the high state, as well 
as the detection of a sharp far-IR spectral break 
in the high state, strongly suggest the existence of 
two synchrotron components: one peaking in the 
mm/sub-mm band and the other peaking in the 
mid-IR band. The component with a peak at the 
mid-IR band can be responsible for emission during 
7-ray flaring states. 

We have applied our leptonic emission model for 
the SEDs during the 7-ray flaring state with a po- 
larization change. Adopting the interpretation of 
the polarization swing involving the propagation of 
the emitting region - that suggested in Paper I - we 
can explain the evolution of the broad-band SEDs 
from Periods D to E during the 7-ray flaring event 
by a shift of the position of the emitting region and 
a change of the viewing angle that are consistent 
with its trajectory. We used the same distribu- 
tion of magnetic fields and only slightly changed 
electron spectra, but required a rather steep strat- 
ification of the external radiation density in the 
form of thermal emission from the dusty torus. 
In this case, the far-IR spectral break requires a 
break in the electron distribution. The observed 
Spitzer-MIPS spectral index 070-160 = 0.35 ± 0.23 
is marginally consistent with the synchrotron emis- 
sion in the fast-cooling regime. 



We also discussed the model in which the 7-ray 
flare is generated within the broad emission line re- 
gion at sub-pc scale from the central black hole ac- 
cording to the jet precession scenario. This model 
explains the mid-IR break during the flaring state 
of Period D by synchrotron self-absorption. Here, 
we also discussed the low-state SED in Period A 
where the mm/sub- mm-band-pcaking synchrotron 
component can be dominant. The model shows the 
related ERC component can explain the steady X- 
ray emission. 
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APPENDIX 

POSSIBLE EXPLANATION OF A LAG BETWEEN THE 7- RAY AND OPTICAL FLARES 

As we discussed in Section |4~T| the multi-band time series imply that during the flaring activity detected in 3C 279 
the optical emission appears delayed with respect to the 7-ray emission. In the context of radiation models adopted 
here (Section [5]), the same electron population produces optical synchrotron photons and also inverse- Compton 7-rays 
in the fast-cooling regime. A lag between the optical and 7-ray flares may therefore result from different profiles of 
decrease of the magnetic energy density u' B (r) and the external (target) radiation energy density u' ext (r) with the 
distance r along the jet, convolved with a non-monotonic profile of the electron injection rate within the outflow. 

In the fast-cooling (FC) regime, the power injected into the relativistic particles P c , inj is immediately radiated away 
and determines the total broad-band luminosity produced by the cooled electrons £ to t,FC- Assuming a strong inverse- 
Compton dominance, i.e., the observed inverse- Compton peak (7-ray) luminosity ~ L 1 being much larger than the 
observed synchrotron peak (optical) luminosity L pt, one has 

Pe.injO") oc Ltot,Fc(r) ^ L 1 {r) + L opt (r) ~ L 1 (r) . (Al) 

while the optical luminosity is 

iopt(r) " L ^ r) K irir 1 ) p +**w < (A2) 

"extV'J a cxtV ) 

where we assumed 5 = Tj being independent on the position r along the jet. Hence, it is clear that while a maximum 
of L^(r) is determined solely by the injection rate P Ci inj(r), a maximum of L opt (r) may in general be quite different, 
depending on particular radial profiles of P e ,lnj(r) and of the ratio [«g/«^ xt ] (r). 

As a specific illustrative example, let us assume that the dissipation region propagating down the jet injects non- 
thermal energy into radiating particles at the rate being a broad Gaussian function of distance r with a maximum at 
ro and a width of ro/v2) 

P Ciinj cx exp [-(r - r ) 2 /rl] , (A3) 
and that magnetic field and external photon field energy densities scale with r as power laws with indices /3b and /3 cxt 

u' B (r) cx r _/?B and Mext( r ) K r~ 0oxt . (A4) 
Then one can find that L 7 has a maximum at r = r , as expected, whereas L opt attains a maximum at 

r CI = y x (1 + Vl + 2(/3 oxt -/3 B )) , (A5) 

which is larger than ro as long as /3 cxt > /3b and thus results in the optical flare lagging the 7-ray flare. This is due 
to the fact that with the magnetic energy density decreasing less rapidly in the jet comoving frame than the external 
radiation energy density, the drop in the injection rate P c , inj(^) between ro and r cr is compensated by the increase in 
the ratio [u' B /u' cxt ] (r). 

Let us further consider the particular values of /3b = 2 and /3 cx t = 4 discussed in Section [5] With such, assuming 
again the electron injection rate being a broad Gaussian function of the distance r along the jet as in the example 
above, the observed time lag between the optical and 7-ray flares At b s can be evaluated as roughly 

It is encouraging that a 10-day lag is expected for Tj ~ 15 and ro — 3pc, which are the bulk Lorentz factor and the 
location of the dissipation region comparable to that inferred from our ERC-IR modeling. 



REFERENCES 

Abazajian, K. N., et al. 2009, 182, 543 

Abdo, A. A., et al. (Fermi-LAT Collaboration) 2009a, 

Astroparticle Physics, 32, 193 
Abdo, A. A., et al. (Fermi-LAT Collaboration) 2009b, ApJ, 699, 

817 

Abdo, A. A., et al. (Fermi-LAT Collaboration) 2010a, Nature, 
463, 919 

Abdo, A. A., et al. (Fermi-LAT Collaboration) 2010b, ApJS, 188, 
405 

Abdo, A. A., et al. (Fermi-LAT Collaboration) 2010c, ApJ, 710, 
810 

Abdo, A. A., et al. (Fermi-LAT Collaboration) 2010d, ApJ, 710, 
1271 

Abdo, A. A., et al. (Fermi-LAT Collaboration) 2010e, ApJ, 721, 
1425 

Abdo, A. A., et al. (Fermi-LAT Collaboration) 2010f, ApJ, 722, 
520 



Abdo, A. A., et al. (Fermi-LAT Collaboration) 2010g, Physical 

Review Letters, 104, 101101 
Abdo, A. A., et al. (Fermi-LAT Collaboration) 2011, ApJ, 726, 43 
Ackermann, M., et al. (Fermi-LAT Collaboration) 2010, ApJ, 

721, 1383 

Ackermann, M., et al. (Fermi-LAT Collaboration) 2012, ApJ, 
751, 159 

Agudo, I., Marscher, A. P., Jorstad, S.G., Larionov, V.M., et al. 

2011a, ApJ, 735, L10 
Agudo, I., et al. 2011b, ApJ, 726, L13 

Albert, J., et al. (MAGIC Collaboration) 2008, Science, 320, 1752 
Aleksic, J., et al. (MAGIC Collaboration) 2011, A&A, 530, A4 
Atwood, W. B., et al. (Fermi-LAT Collaboration) 2009, ApJ, 697, 
1071 

Baars, J. W. M., et al. 1977, A&A, 61, 99 
Begelman, M. C, & Sikora, M. 1987, ApJ, 322, 650 
Bessell, M. S., Castelli, F., & Plez, B. 1998, A&A, 333, 231 
Blandford, R. D., & Levinson, A. 1995, ApJ, 441, 79 
Bloom, S. D., & Marscher, A. P. 1996, ApJ, 461, 657 



Radio to high-energy 7-ray observations of 3C 279 in 2008-2010 



23 



Bock, D.C.-J., ct al. 2006, Proc. SPIE, 6267, 13 
Bonning, E. W., et al. 2009, ApJ, 697, L81 
Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 
245 

Chatterjee, R., et al. 2008, ApJ, 689, 79 

Collmar, W., et al. 2001, Gamma 2001: Gamma-Ray 

Astrophysics, 587, 271 
Collmar, W., Bottcher, M., Krichbaum, T. P., et al. 2010, A&A, 

522, A66 

Cutri, R. M., ct al. 2003, VizieR Online Data Catalog, 2246, 
Daly, R.A., & Marscher, A.P. 1988, ApJ, 334, 539 
D'Ammando, F., et al. 2011, A&A, 529, A145 
Dermer, C. D., Schlickeiser, R., & Mastichiadis, A. 1992, A&A, 
256, L27 

Dermer, C. D., & Schlickeiser, R. 1993, ApJ, 416, 458 

Eachus, L. J., & Liller, W. 1975, ApJ, 200, L61 

Giuliani, A., et al. 2009, A&A, 494, 509 

Greiner, J., et al. 2008, PASP, 120, 405 

Gu, M., Cao, X., & Jiang, D. R. 2001.MNRAS, 327, 1111 

Edelson, R. A., & Krolik, J. H. 1988, ApJ, 333, 646 

Edelson R., ct al. 2002, ApJ, 568, 610 

Engelbracht, C. W., et al. 2007, PASP, 119, 994 

Gordon, C. W., et al. 2007, PASP, 119, 1019 

Hartman, R. C, et al. 1992, ApJ, 385, LI 

Hartman, R. C, et al. 1999, ApJS, 123, 79 

Hartman, R. C, et al. 2001a, ApJ, 553, 683 

Hartman, R. C, et al. 2001b, ApJ, 558, 583 

Hermsen, W., et al. 1993, A&AS, 97, 97 

Homan, D. C, Lister, M. L., Kellermann, K. L, et al. 2003, ApJ, 
589, L9 

Jorstad, S. G., et al. 2004, AJ, 127, 3115 
Jorstad, S. G., et al. 2005, AJ, 130, 1418 
Jorstad, S. G., et al. 2010, ApJ, 715, 362 

Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, 
A&A, 440, 775 

Komissarov, S.S., & Falle, S.A.E.G. 1997, MNRAS, 288, 833 
Koyama, K., Tsunemi, H., Dotani, T., et al. 2007, PASJ, 59, 23 
Kriihler, T., ct al. 2008, ApJ, 685, 376 
Larionov, V. M., et al. 2008, A&A, 492, 389 
Lynds, C. R., Stockton, A. N., & Livingston, W. C. 1965, ApJ, 
142, 1667 

Mannheim, K., & Biermann, P. L. 1992, A&A, 253, L21 
Maraschi, L., Ghisellini, G. & Celotti, A. 1992, ApJ, 397, L5 
Maraschi, L., ct al. 1994, ApJ, 435, L91 
Mattox, J. R., et al. 1996, ApJ, 461, 396 
McNaron-Brown, K., et al. 1995, ApJ, 451, 575 
Mitsuda, K., et al. 2007, PASJ, 59, 1 



Moderski, R., Sikora, M., & Blazcjowski, M. 2003, A&A, 406, 855 

Mor, R., & Netzer, H. 2012, MNRAS, 420, 526 

Morrison, R., & McCammon, D. 1983, ApJ, 270, 119 

Nalewajko, K., 2010, IJMPD, 19, 701 

Nalewajko, K., & Sikora, M. 2009, MNRAS, 392, 1205 

Nandra, K., ct al. 1997, ApJ, 476, 70 

Nandikotkur, G., et al. 2007, ApJ, 657, 706 

Nilsson, K., Pursimo, T., Villforth, C, Lindfors, E., & Takalo, 

L. O. 2009, A&A, 505, 601 
Ogle, P. M., Wehrle, A. E., Balonek, T., & Gurwell, M. A. 2011, 

ApJS, 195, 19 
Pian, E., et al. 1999, ApJ, 521, 112 
Poole, T. S., et al. 2008, MNRAS, 383, 627 
Poutanen, J., & Stern, B. 2010, ApJ, 717, L118 
Raiteri, C. M., Villata, M., Lantcri, L., Cavallonc, M., & Sobrito, 

G. 1998, A&AS, 130, 495 
Raiteri, C. M., et al. 2011, A&A, 534, A87 
Reach, W. T., et al. 2005, PASP, 117, 978 
Richards, J. L., et al. 2011, ApJS, 194, 29 

Sault, R.J., Teuben, P.J., & Wright, M.C.H. 1995, in ASP Conf. 
Ser. 77, Astronomical Data Analysis Software and Systems IV, 
cd. R.A. Shaw, H.E. Payne, & J.J.E. Hayes (San Francisco: 
ASP), 433 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 
525 

Sikora, M., Bcgclman, M. C, & Rees, M. J. 1994, ApJ, 421, 153 
Sikora, M., Stawarz, L., Moderski, R., Nalewajko, K., & 

Madejski, G. M. 2009, ApJ, 704, 38 
Sikora, M. 2011, IAU Symposium, 275, 59 

Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 
1163 

Stansberry, J. A., et al. 2007, PASP,119, 1038 
Takahashi, T., et al. 2007, PASJ, 59, 35 

Tanaka, Y. T., Stawarz, L., Thompson, D. J., et al. 2011, ApJ, 
733, 19 

Tody, D. 1993, Astronomical Society of the Pacific Conference 
Series, 52, 173, editor: R. J. Hanisch, R. J. V. Brisscndcn, & 
J. Barnes 

Ulrich, M.-H., Maraschi, L., & Urry, C. M. 1997, ARA&A, 35, 445 
Vaughan, S., et al. 2003, MNRAS, 345, 1271 
Villata, M., ct al. 1997, A&AS, 121, 119 
Villata, M., ct al. 2008, A&A, 481, L79 
Villata, M., ct al. 2009, A&A, 504, L9 

Watanabe, M., et al. 2005, Pub. Astron. Soc. Pacif., 117, 870 
Wehrle, A. E., et al. 1998, ApJ, 497, 178 

Wolff, M. J., Nordsieck, K. H., & Nook, M. A. 1996, AJ, 111, 856 
Woo, J.-H., & Urry, C. M. 2002, ApJ, 579, 530 



